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POLYNOMIALS  WITH  RESTRICTED  COEFFICIENTS  AND  THEIR  APPLICATIONS 


PROMETHEUS  INC. 


ABSTRACT 

Certain  design  restrictions  growing  out  of  antenna  theory 
yield  a  beautiful  class  of  complex  variables  problems.  Our  work  was 
devoted  to  formulating  these  problems  in  mathematical  terms,  solving 
some  of  them,  and  beginning  work  on  the  others.  One  important  result 
achieved  was  the  development  of  a  new  method  of  estimating  6aussian- 
type  exponential  sums.  Improvements  of  our  previous  results  in  null 
steering  and  notch  filtering  were  also  attained.  Other  findings  were 
obtained  in  the  areas  of  5  the  robustness  of  polynomials  with 
unimodular  coefficients;  the  effect  of  errors  in  such  standard  and 
crucial  approximations  as  the  far-field,  Fresnel,  and  Doppler 
compensation;  the  effect  of  errors  in  the  noise  covariance  matrix;  and 
the  Parabolic  Equation  Method  in  underwater  acoustics. 
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STATEMENT  OF  WORK 


1.  Convert  the  Kahane  "Gauss  Sums  +  Probabilistic  Choices"  proof  into 
a  constructive  one,  and  apply  these  polynomials  to  the  design  of  low 
peak-factor  signals  and  other  engineering  problems. 

2.  Investigate  other  possible  methods  of  constructing  Kahane-type 
polynomials. 

3.  Investigate  employing  the  Byrnes  and  Kahane  polynomials  to  design 
new  reflection  phase  gratings,  and  consider  the  application  of  these 
ideas  to  the  quieting  of  an  objects  response  to  radar  and  active 
sonar,  and  the  design  of  baffles  used  to  quiet  machinery  noise  from 
submarines. 

4.  Consider  the  incorporation  of  the  Prometheus  methods  of 
analytically  choosing  shading  coefficients  into  adaptive  antenna 
problems. 

5.  Determine  the  applicability  of  Kahane’s  methods  to  the  Erdos 
problem  involving  polynomials  with  coefficients  +1,  and  investigate 
the  minimum  L*  norm  of  such  polynomials. 

6.  Investigate  the  minimization  of  the  peak-factor  in  an  inverse 
Fourier  transform,  given  the  amplitude  spectrum. 

7.  Conduct  research  into  whether  polynomials  based  upon  the  Byrnes 
construction  can  be  used  to  produce  notched  filters  with  more  than  one 


notch. 


I.  INTRODUCTION 

Polynomials  with  restricted  coefficients  have  been  objects  of 
intense  interest  throughout  twentieth  century  mathematics  and 
engineering.  Despite  the  large  amount  of  effort  which  has  been 
directed  towards  analyzing  their  properties,  many  deep  and  fascinating 
questions  remain.  We  continue  our  attack  upon  these  problems, 
employing  both  traditional  and  novel  methods.  This  report,  with  its 
appendices,  presents  the  details  of  our  recent  work. 

We  have  seen  that  certain  design  restrictions  growing  out  of 
antenna  theory  yield  a  beautiful  class  of  complex  variables  problems. 
Moreover  these  problems  form  an  integral  part  of  classical 
mathematical  analysis,  as  the  study  of  such  polynomials  was  initiated 
by  G.  H.  Hardy  L46,p.  1993,  and  furthered  by  J.  E.  Littlewood  C  28-30  3, 
P.  Erdos  C203,  and  many  others.  Our  research  continues  to  focus  upon 
both  aspects  of  this  remarkable  intertwining  between  the  disciplines 
of  pure  mathematics  and  engineering. 

The  great  success  of  J.  P.  Kahane  C 26 3  in  solving  the 
Livtlewcod  conjecture,  by  showing  that  there  indeed  exist  polynomials 
with  unimodular  coefficients  whose  modulus  is  essentially  constant  on 
the  unit  circle,  lead  us  to  reconsider  the  Erdos  variant  of  this 
problem.  Namely,  do  such  polynomials  exist  when  the  coefficients  are 
further  restricted  to  be  +  1?  Kahane’s  answer  to  Littlewood’s  question 
completed  a  series  of  researches  by  Seller  C2-53,  Byrnes  C  13,  143, 
Kort.er  C  27  3,  Newman  C3-5,  34  3,  and  Littlewood  himself.  As  just 
mentioned  the  answer  was  an  unqualified  yes.  near  constancy  for 
polynomials  with  unimodular  coefficients  is  definitely  available.  In 
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contrast  to  this  there  has  been  very  little  progress  in  answering 
Erdos’  question,  and  this  has  been  one  focus  of  our  research. 

We  are  conducting  a  two  pronged  attack  upon  the  Erdos 
question.  One  tack  we  are  taking  is  to  attempt  to  exploit  the 
breakthrough  of  Kahane,  which  we  believe  was  due  to  his  ingenious  use 
of  randomness  and  probability  in  his  construction.  Behind  his  and 
previous  approaches  was  the  idea  of  Gauss,  viz.  the  “Gauss  Sums.*  To 
put  it  quite  simply  we  feel  that  Littlewood’s  problem  was  vanquished 
by  the  “equation* 

Kahane=Gauss  Sums  +  Probabi 1 istic  Choices 
Since  the  Gauss  Sums  can  also  be  written  as  polynomials  with 
coefficients  +1,  we  devoted  considerable  study  to  such  sums.  Our 
entirely  new  approach  to  estimating  them,  together  with  the  results 
achieved  thusfar,  are  described  below. 

The  above  describes,  to  some  extent,  the  first  prong  of  our 
attack  on  the  Erdos  problem.  The  second  is  our  approach  to  the 
possibility  of  a  negative  answer  to  the  Erdos  question.  Indeed 
numerical  evidence  seems  to  favor  this  choice.  A  computational  scheme 
based  upon  the  Bose-Einstein  model  for  Statistical  Mechanics  has 
produced  calculations  of  autocorrelations  which  seem  to  insist  that 
the  L*  norm  of  a  +  1  type  polynomial  is  >  (6/5  > 1 '■*<  n  >  * /:a,  and  this  if 
true  would  surely  yield  the  ‘no"  answer  to  the  Erdos  problem.  Since 
this  Bose-Einstein  model  can  be  followed  up  theoretically,  and  not 
only  numerically,  we  envision  a  possible  proof  of  the  negative  answer, 
if  such  is  the  case.  What  we  have  proven  thusfar  is  that  the  norm 
of  the  well-known  Shapiro  polynomials  of  type  is  ( 4/3 )*'"*( n )  1 
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asymptotically.  Details  of  these  results  appear  below,  and  in 
Appendix  B. 

The  locations  of  zeroes  of  polynomials  with  restricted 
coefficients  is  of  interest  from  both  a  theoretical  and  a  practical 
point  of  view.  Under  some  circumstances,  such  as  those  encountered  in 
peak  power  limited  transmitting  C4Q1,  the  design  of  low  peak-factor 
C393  and  low  crest  factor  C73  signals,  and  the  design  of  digital  test 
signals  C361,  one  wants  the  modulus  of  the  polynomial  to  be  as  close 
to  constant  as  possible  on  the  unit  circle.  It  seems  reasonable  that 
such  a  condition  would  be  achieved  by  polynomials  whose  zeroes  are  as 
far  as  possible  from  the  unit  circle.  Certain  applications,  on  the 
other  hand,  require  zeroes  at  specific  locations  on  the  unit  circle. 
For  example,  the  classical  mathematical  problem  in  notch  filter  design 
is  to  produce  a  polynomial  whose  magnitude  on  the  unit  circle  is  close 
to  constant  in  almost  all  directions,  but  which  has  a  small  number  of 
deep  nulls  < “notches")  at  specified  points.  In  C163  the  construction 
of  C133  is  employed  to  produce  such  a  polynomial  having  one  null,  with 
the  added  feature  that  all  coefficients  have  the  same  magnitude.  We 
have  expanded  upon  these  methods  to  produce  multiple  notches,  while 
maintaining  the  near  constancy  of  the  magnitude  of  the  polynomial  at 
all  points  on  the  unit  circle  excluding  small  neighborhoods  of  the 
notches.  Again,  the  results  appear  below. 

Mull  steering  is  a  crucial  application  which  requires  locating 
zeroes  at  specific  points  on  the  unit  circle.  There  are  several 
factors  which  must  be  considered  in  the  design  of  null  steering 
algorithms.  In  addition  to  the  basic  problem  Of  placing  the  nulls  the 
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main  beam  must  be  steered,  the  width  of  the  main  lobe  controlled,  and 
the  sidelobe  levels  must  be  sufficiently  below  that  of  the  main  lobe. 
Control  of  the  sidelobe  level  is  usually  achieved  by  attenuating  the 
shading  coefficients  as  one  moves  away  from  the  center  of  the  array. 
Often  these  attenuation  factors  (Chebyshev,  Taylor,  etc.)  are  chosen 
in  advance,  and  may  not  be  easily  altered  once  the  array  is  in  place. 
This  leads  directly  to  a  beautiful  mathematical  question,  similar  to 
the  peak-factor  problem  in  engineering  discussed  earlier,  which  we 
have  attacked  by  the  methods  employed  to  study  polynomials  with 
coefficients  of  magnitude  one: 

Given  the  magnitude  of  the  coefficients  of  a  polynomial  P,  a 
finite  subset  S  of  the  unit  circle  C,  and  a  point  p€C  distinct 
from  those  in  S,  choose  the  phases  of  these  coefficients  so 
that  P(z)^0  for  all  2<S,  the  maximum  on  C  of  iP(z)l  occurs  at 
zap,  and  the  maximum  of  1 P <  z ) I  on  a  subset  of  C  excluding  an 
appropriate  interval  < the  beamwidth)  around  p  is  as  small  as 
possible. 

Currently  the  most  widely  used  class  of  null  steering  methods 
is  known  as  adaptive  nulling  C 1 , 9, 24, 25, 32, 45 3.  Adaptive  arrays  have 
developed  over  the  past  twenty-five  years  as  the  preferred  method  of 
reducing  the  performance  deter iorat ion  in  signal  reception  systems 
which  is  inevitably  caused  by  undesired  noise  entering  the  system. 
Sources  for  this  noise  include  multipath  affects,  electronic 
countermeasures,  clutter  scatterer  returns,  antenna  location  errors, 
array  ele«*ent  thermal  noise,  etc.  The  prol  i  f erat  ion  of  such  noise 
sources  has  greatly  increased  the  importance  of  interference 


suppression  in  essentially  all  applications.  Although  such  adaptive 
methods  as  the  Widrow  least  mean  squares  (LMS)  and  Howel ls-Applebaum 
sidelobe  canceller  have  achieved  considerable  success,  difficult 
problems  remain.  Foremost  among  these  are  poor  transient  response, 
signal  cancellation  resulting  from  interaction  between  signal  and 
interference,  excessive  computation  time,  and  sidelobe  degradation 
uihen  jammer  cancellation  is  attempted.  A  secondary  problem  is  the 
lack  of  control  in  adaptive  algorithms  of  the  dynamic  range  of  the 
weights. 

These  methods  are  indirect  adaptive  schemes,  they  do  not 
explicitly  form  an  estimate  of  the  directions  of  arrival  of 
interfering  sources  or  explicitly  steer  nulls  in  those  directions.  A 
scheme  in  which  these  two  tasks  are  actually  performed  can  be  called  a 
direct  adaptive  algorithm.  Thus  one  approach  to  the  solution  of  such 
problems  is  to  complement  an  appropriate  indirect  adaptive  algorithm 
with  the  analytic  null  steering  methods  described  in  1173.  In  this 
way  the  actual  noise  suppression  achieved  can  be  enhanced  beyond  that 
which  would  be  available  through  either  adaptive  or  analytic  methods 
exclusively.  We  have  begun  to  explore  the  possibility  of  this  'Direct 
Adaptive  Antenna  System,*  and  will  continue  this  work  during  Phase  II, 

Other  questions  which  have  been  the  focus  of  our  research 
effort,  and  which  we  report  on  below,  include:  the  robustness  of 
polynomials  with  unimodular  coefficients;  the  effect  of  errors  in  such 
standard  and  crucial  approximations  as  the  far-field,  Fresnel,  and 
Doppler  compensation;  the  effect  of  errors  in  the  noise  covariance 
matrix;  and  the  Parabolic  Equation  Method  in  underwater  acoustics. 


I I .  WEYL  SUMS 


In  connection  with  the  beamforming  problems  under 
consideration,  uie  were  required  to  obtain  certain  trigonometric  sum 
estimates.  In  particular,  sums  such  as 

•> 

g  i«M- 

k«l 

made  their  appearance,  and  we  desired  " Gauss i an- 1 ike"  estimates  for 
these  "Gaussian-1 ike"  sums.  In  short,  we  wanted  the  bound  of  0  vTrT  for 
said  sums  taken  over  intervals,  I,  of  length  of  order  n  and  for  t  of 
the  order  of  1/n. 

The  available  method  for  handling  such  sums  is  that  of  Hermann 
Weyl.  Realizing  the  intractability  of  sums  of  exponentials  of 
quadratics,  he  reduced  such  sums  to  exponentials  of  1 inears  (geometric 
sums!)  by  simply  multiplying  by  their  conjugates.  Unfortunately, 
however,  his  method  leads  to  bounds  like  0  Jn ~logn  instead  of  the  0  JrT 
that  we  desire. 

To  see  the  emergence  of  this  logn  factor,  let  us  examine  his 
procedure  in  some  detail.  Assume  that  I  is  an  interval  of  length 
and  write 


Z2 

e ***•  so  that 

kd 

|S|»  »  SS  =  5^ 

Ml 


2  V-  2  c —  2  2 
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If  we  then  write  k-j=v,  so  that  K+j=2k-vf  and  call  l._  the  intersection 
of  the  two  intervals,  I  and  I  translated  by  <-v,  we  obtain 


*VI  2 

(1.)  |s|a=£_  e 

V--M  keXv 


3  A  tvk 


But  here  is  the  rub:  one  has  the  choice  of  throwing  away  the 
factor  e-1*"'  by  “clashing  through”  and  thereby  losing  accuracy,  or 
being  stuck  with  the  same  intractable  sum  of  exponentials  of 
quadratics.  Weyl  opted  tor  the  ‘crashing  through."  The  resulting  bound 
is  (with  n^=  1 1 vj  ) , 


n-l 

H 


n-l 


0  \ 

sin  tv  -  I sig ti// 


I 


■/zi 


~rr  n'f  t r 

<  n+— --V-  (if  o<t<- — ■,  e.g.) 

-£  L _  V  2-h 

1/-I 


and  thereby  emerges  the  onerous  logn. 

Our  approach,  then,  is  to  make  the  other  choice  and  face  the 
troublesome  ,  the  trick  being  that  this  sum,  thaugt  still 

intractable,  can  be  directly  related  to  our  original  sum.  An  estimate 
of  the  desired  quantity  is  given  ig  terms  of  itself,  and  this  will 
prove  effective. 

So  let  M  be  the  exact  bound  for 

keX 


over  all  intervals,  I,  of  length  <n  and  all  real  quadratic 
polynomials,  P(k>,  with  leading  coefficient  t#o.  We  will  first  treat 
the  restricted  case  where  o<t<TT/2n.  The  general  case  will  follow  from 
this  one  since  this  says  n^  ~TF /2t  so  that  a  general  sum  obviously 
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splits  into  pnt/'ffj  of  these  restricted  ones.  Our  aim  is  to  shoui  here 
that  M=0( t-1'2)  and  this  is  the  promised  'Gauss-like"  bound  since  our 
real  interest  is  in  the  case  of  t  of  exact  order  1/n. 

To  begin  with  we  need  the  following  well  known  facts. 

Lemma  1  (Summation  by  parts): 


n 


akbk  <<  PS(a,<)  V(bfc)  where  PS  (partial  sum  bound)  is 

bl 


Max  \T  a* 

htn  IfeTnx 


,  and  V  (total  variation)  is 
0-1 


min(  |b 4,  |t>4>  +  y~  j  . 

k-\ 

(In  particular,  if  the  bw  are  monotomc  and  of  the  same  sign,  then 
0<b*>  =  Max  (  |  bi|,  |b„|  ) ) . 

Lemma  2  (Gibbs’  phenomenon)  The  partial  sums  of  the  series 


oo 

5(h  k-o- 


are  uniformly  bounded  by  2  (not  by  Tf/2,  which  might  be  guessed  since 
this  is  the  Fourier  series  of  <TF-0>/2  on  o<0<,2 IT  ). 

Letting  I=C  r  +  1 ,  r+ro  3,  where  m^n,  the  same  algebraic 

manipulations  as  before  lead  to 


Z 


n-l 

*m-*- 1  m  J 
*1 


-(1 t^-cQv  f[v 

_ _ 

*  in 


) 
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with  B  —  a+(2r+l)t,  A  =  a+(2r+2m+l ) t . 


To  estimate  the  right  hand  sum,  we  split  it  into 


i /<-f 


and  S3 


Jit  < 


For  S2,  we  use  our  Lemma  1,  with  the  av  the  numerators,  and  the  bv= 
1/sin  vt.  The  PS  term  is  clearly  bounded  by  2M  and  the  V  term,  by 
monotonicity,  is  bounded  by 


Altogether,  then,  we  get 

Sa  <<  "lyMr-1"* 

To  bound  St,  u  different  tack  must  be  taken.  We  break  S, 
further  into  the  two  sums 

Im  S  ^ — — . -  and  Im*V  — _ _ _ _ 

fa*  Vt  **  Vt 

In  both  cases,  the  dropping  of  the  quadratic  term  in  the  exponent 
causes  an  error  or  <<vat  in  each  term  of  the  numerator.  Thus,  each  of 
these  sums  is  replaceable  by  a  sum 


T  *  T  ~~r~ 

* —  nn  h 

✓<f- 1 


with  an  error 
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Finally,  we  must  effect  the  estimation  of  our  two  sums  T. 
This  is  achieved  by  writing  it  as 


T  =4- 


Srh  VI 


*  v 


/  -Sin  vC  (— - r 


The  first  of  these  is  bounded  by  2/t  by  Lemma  2,  and  crude  methods 
suffice  for  the  second  sum.  Since 


sin  x  % 


<  x  (or  sin  x  >• 


in  the  interval  (o,17"/2)  (fun  for  the  reader?),  we  obtain,  by  crashing 
through,  the  estimate 


vt<l  for  this  second  sum. 


Altogether,  then,  these  many  bounds  combine  to  give 


NT-  X  |*“  “tT^TT 

j  <  n+*7T  Mt  +  +2, 


k*T 


and  by  picking  a  and  m  Judiciously,  this  means  that 


9  IT 

Ma<n-»  M. 

t  yft 


Hence,  as  promised,  we  have  estimated  our  bound,  M,  in  terms 


of  itself.  Finally,  then,  we  write  this  as 
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(M-  -^=)  <  n+  and  obtain  M  <  r  +  [ n+^~  <~Tt 
•jJT  t  ~  v  t  ~  'Jt 

which  is  what  we  desired  for  the  restricted  sum.  For  the  general  sum 
we  need  only  multiply  by  and  obtain 

M<  +4n\Tt  - 

Indeed,  we  have  saved  our  precious  logarithm! 

Extension  of  the  above  results  to  the  case  where  P(k)  is  a 
third  degree  polynomial,  and  then  to  arbitrary  polynomials,  will  be 
pursued  during  Phase  II. 
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III.  L*  NORMS 


Investigations  were  continued  into  the  L*  norm  of  polynomials 
with  coefficients  +1.  Revisions  in  the  previous  paper  C3S3  were  made, 
and  the  revised  paper,  which  is  Appendix  B  in  this  report,  will  appear 
in  the  American  Mathematical  Monthly. 

To  describe  the  new  results,  we  again  employ  the  notation  of 

Appendix  B.  Thus, 

o-i 

P(z)  =  each  <»,=1  or  -1,  z=ears’1®‘,  0<-*Kl. 

o 

As  shown  previously, 

.  q 

j/PlJ  “  2_.  <».  <m- 

l’ 4  j  ♦ X; 

A  straightforward  computation  yields  the  lower  bound  na+2Cn/2  3  for 
this  sum.  Also,  by  considering  -P(z)  and  P(-z),  it  is  clear  that, 
without  loss  of  generality,  it  may  be  assumed  that  Using 

this,  an  exhaustive  computer  search  for  the  minimum  L*  norm  was  made 
for  n<22,  on  an  IBM-XT  compatible,  with  the  results  appearing  in  Table 

1.  Note  that  the  lower  bound  na+2Cn/23  is  achieved,  with  a  unique 

combination  of  <’s,  for  n=3,5,7, 1 1, 13,  and  it  is  also  achieved,  but 
not  uniquely,  for  n=4.  Also  observe  that  the  number  of  computations 
required  to  evaluate  the  sum  grows  very  rapidly  with  n,  and  the 

exhaustive  search  becomes  impossible,  even  for  the  largest  current 

computers,  when  n  reaches  about  30.  Thus,  other  methods  are  necessary 
to  decide  the  minimum  L*  norm  of  such  polynomials  and  to  possibly 
answer  the  Erdos  question  mentioned  in  C353.  Work  on  this  problem  is 
ontinuing. 
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IV.  NOTCH  FILTERS  AND  NULL  STEERING 


Investigations  were  also  continued  into  notch  filters 
employing  coefficients  of  equal  magnitude.  Revisions  in  the  previous 
paper  C163  were  made,  and  the  revised  paper,  which  is  Appendix  D  in 
this  report,  has  been  tentatively  accepted  (pending  approval  of  the 
changes)  for  publication  in  the  IEEE  Transactions  on  Acoustics,  Speech 
and  Signal  Processing. 

The  new  results  in  this  area  involve  the  design  of  notch 
filters  with  multiple  notches,  but  still  with  dynamic  range  equal  to 
1.  The  basic  idea  is  to  again  begin  with  the  function  P(ff)  of  [1  1, 
but  instead  of  removing  the  first  N  terms  (i.e.,  all  terms  with  k=0), 
as  was  done  in  C 1  3  to  obtain  all  terms  for  several  values  of  k 
are  removed.  This  results  each  time  in  a  notch  filter,  with  the 
number  of  notches  equal  to  the  number  of  k’s  that  are  removed.  Also, 
as  all  nonzero  coefficients  have  magnitude  one,  these  remain  filters 
with  unit  dynamic  range.  The  worth  of  this  method  may  be  seen  in  the 
graphs  on  the  next  18  pages,  which  exhibit  filters  of  degree  Q  to 
2824,  with  from  1  to  4  notches. 

Another  previous  paper  £  1 73,  dealing  with  the  important 
subject  of  null  steering,  required  revision.  This  was  accomplished  as 
part  of  the  current  effort,  and  the  resulting  paper,  which  is  Appendix 
E  in  this  report,  will  appear  in  the  February  1988  issue  of  the  IEEE 
Transactions  on  Antennas  and  Propagation. 
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A  Single  Notch  Filter,  Degree  224 
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A  Double  Notch  Filter,  Degree  483 
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V.  APPROXIMATION  ERRORS 


Another  problem  considered  was  the  development  of  new  methods 
of  incorporating  Doppler  into  standard  models  of  a  target  moving  with 
respect  to  an  array.  In  order  to  make  the  question  mathematically 
tractable,  the  simple  case  of  a  target  path  parallel  to  a  passive  line 
array  of  equally  spaced  elements  was  analysed.  Generalizations  of  the 
concepts  described  herein  will  be  explored  in  Phase  II. 

Thus,  suppose  a  moving  target  is  generating  a  signal  s(t) 
which  is  being  received  by  an  N-element  line  array  with  interelement 
spacing  d.  Let  r^,(t>  be  the  distance  at  time  t  between  the  target  and 
the  n-th  array  element  A„(l<n<N),  y  be  the  perpendicular  distance 
between  the  target  path  and  the  line  of  the  array,  c  be  the  speed  of 
the  signal,  v  be  the  speed  of  the  target,  and  T„(t)  -  r„(t)/c  be  the 
time  delay  from  the  target  to  A*,.  Neglecting  Doppler  (and  noise)  for 
the  moment,  the  received  signal  at  A„  is  a„ ( t  )=s (t-T„<  t  !)/r„(  t ) . 

The  Doppler  effect  can  be  represented  by  a  multiplicative 
factor  u  in  the  time  variable,  so  that  s  above  becomes  s(ut-Tn(ut )) . 
If  u  were  constant,  the  required  transformation  of  this  factor  into 
frequency  domain  would  be  simple,  namely,  the  Fourier  transform  S(w) 
of  stt)  would  become  <  l/u.'S(w/u) .  However,  in  actuality,  u  is  a 
function  of  both  n  and  t,  u=l+frt/c.  Here,  r„  represents  the  rate  of 
change  of  r„(t)  with  respect  to  time  and  is  negative  when  the  target 
is  approaching  A„.  One  approximate  yet  reasonable  method  of  dealing 
with  the  variable  Doppler  ir.  to  assume  that  for  each  n,  u  changes 
exactly  once  as  the  target  moves,  with  this  change  occurting  at  the 
time  tQ  when  the  tarqat  is  exactly  opposite  the  middle  array  element. 
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Letting  x(t)  be  the  target  location  at  time  t  and  n * =C (N+l ) /2 3,  the 
situation  is  given  in  figure  1.  „  A 


From  this  figure,  clearly 
x<t>=x<0)+vt  and 
r„*  ( t )»y*+(x ( t )— x ( t„)-<n-n*  )d)*. 

To  obtain  a  potentially  useful  estimate  of  r„,  so  as  to 
estimate  u,  let  c T *— 1  if  t<t»,  f£~ =+l  if  t> t„,  be  the  length  of  the 
time  interval  over  which  the  signal  is  measured,  and  r  be  the  average 
value  of  all  rn(t), 

,  ni  ,  r* 

r„( t )dt. 


Applying  the  Mean  Value  Theorem  to  the  above  equation  for  rn(t)  yields 


where  r»  =  y*+  ( x ( t >-x < t0 n*  ) d)*  and  is  between  n  and  n*. 

Approximating  ^  by  <n+nM/2  and  using 

r„-r"w=  < n-n *  )d^(n-n*  )d-2(x<t>-x(t„))^  gives 

\  * 

Y'  -  V  ,  '  t,  C r -  rt  .  A  ~  ^  ^9  W  ^  -^r—Cn—n'i^- 

v 

Since  r£.  <t>=ya+(x<t)-x<t.))*. 


Combining  the  above,  we  see  that 


Finally,  letting  u~  be  our  approximation  for  u,  so  that 
u-*H~rCcf  /l- 4_+  4-  (n‘-n)dl, 

C  1/  X  Z  T 

v  r  r 

the  Fourier  transform  of  s(t)  now  becomes  ( l/u*)SCw/u*  ),  and  this 
quantity  is  employed  in  the  processing  of  the  received  signal. 

New  methods  of  accounting  for  the  effect  of  wavefront 
curvature  upon  beamforming  were  also  considered.  A  preliminary  result 
is  given  below,  and  this  effort  will  continue  in  Phase  II. 

Referring  to  figure  the  receiving  array  is  again  a  line  array 
of  N  equally  spaced  elements,  with  interelement  spacing  d.  The  distance 
from  a  particular  source  to  the  n-th  receiving  element  is  denoted  by  r+x,, 

(with  x,*0),  and  f  is  the  indicated  angle  from  the  first  array  element  to 
the  source. 


SwCt, 


Figure  3k 
-3<»- 


(r+:<B  )*=r*  +  (  (n-1 )  d)*  -2r  (n-1 )  d  cos  p. 


(1) 


Assuming  that  the  array  length<<r  and-  neglecting  all  terms  beyond  the 
first  order,  (1)  yields 

;<„*(  (n-1 )  d)*  /2r  -  (n-l)d  cos  p.  C  2.  ) 

Letting  F^,  be  the  voltage  response  of  the  n-th  element,  and  normalizing 
so  that  RT»lt  it  follows  that  R„ =r/ Cr+x„ ) .  Now,  with  v*  denoting  the 
output  voltage  and  a„  the  shading  coefficient  of  the  n-th  element,  and 
assuming  a  sinusoidal  input  of  frequency  w,  it  follows  that  the  output 
voltage  of  the  array  is 

V  U 

V-5""a*v*®/  a,y  <r/(r+x*  ) )  exp(jwt)  exp(2TT  jx*  /X  / .  <3  ) 

n»T 

Here  A  is  the  wavelength  of  the  source,  so  that  2TTx* /A  is  the  phase 
delay  from  the  first  to  the  n-th  array  element. 

In  most  cases  values  for  V  will  be  obtained  computationally,  as  a 
closed  form  for  the  sum  in  c "S  )  is  usually  impossible  to  obtain. 
However,  if  the  shading  coefficients  a*  are  all  1,  and  if  it  is  assumed 
that  the  x„  increase  linearly  from  x,  to  )^,  then  the  closed  form 
solution  for  V  is 


where  ^|=*(N-1 )  d* /2r  -  d  cos  p. 

To  derive  (A),  observe  that  the  assumption  that  increases 
linearly  from  x,  (=0>  to  x„ ,  when  combined  with  <2),  yields 
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*->  --  jrr  xn  *  M  (-57^-  ■  M  y„ , 

with  defined  as  in  .  From  >,  the  fact  that  all  a„=l, 

assumption  that  array  length<<r,  follows 

j^?Myv  •  x7r 


(.5-) 


»nd  the 


V  '  — 

h  r* 


i rH/»  ,w< 


cw 


Defining 


1 


«•£ 


iL  jury,/*-'] 


i-e,  Jutyh* 


n=l 


( 1> 


it  follows  from  (b)  and  the  termwise  differentiation  of  (7)  that 


Since 


--  £  v„  *  eJ^  (3  (f)+  U-  jYtjJ . 

n-/ 


o'M- 

•J  I  -  * 


i« 


V 


5^ 


(wi 


(&> 


(.4-  •>  is  an  immediate  consequence  of  (b)  ,  (l)  ,  and  (fi  )  + 


Figure 'J*  Configuration  for  the  Fresnel  Approximation 
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■V  ni  rf  /<  j*  vy>  Vh  r> 


Next  we  analyze  the  error  in  the  Fresnel  Approximation 


which,  for  R-r  in  figure  3,  is 

R-r?sv  sin0*Kvacosae-) /2r. 

This  approximation,  useful  when  v<<r  (and  even  more  useful  when  •O-is 
small),  arises  by  truncating  the  power  series  expansion  of  the  square 
root  in  the  equation 


R=r*y  1+pc  +-~r- 


sin 


To  estimate  the  error  in  this  approximation,  we  require  more  detail  in 
the  diagram,  as  in  figure  4. 


By  the  law  of  sines, 

v/sin<  ftt+'i') /2)  =  ( R-r  )/sin ,  or 
v*/cos<'t'/2)=!(R-r >/(sinO“  cosrlV2+casO-  sin't//2) ,  so 
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R-r  =  v  sin©- +  v  cosO*  tan'i'/2=v  sirr©-  +  v  cos©~  ( 1-coS^O /siriir'.  (*) 

From  the  law  of  cosines,  cos  'ir'  =(ra+R2-va)  /2rR,  and  from  the 
law  of  sines,  v/sinT=R/sin <-©+H72)=R/cos©;  so  sirti^vcosO/R.  Inserting 
these  in  the  expression  for  R-r,  we  get 


Denoting  the  desired  error  by  E,  we  have 

E=R— r-v  5in©wv2co5:W2r 

y-C.R-r^ 

=R-r-v  sin©'  -  Z^Z “  +  — — -vacos*©/2 r 


°va  ( 1— cos*©**  <  R-r )  a/va) /2r 

v  fts-  r  .  ft  -  v- 

=  -  (sin© — )(sin©-  +  ~~ — ) 

XT'  V  v 

=  — 7f  (-cos©-  tanv4V2>  <2sin©*  +  cos©"  tarfi'/2, 

>v- 

where  the  last  equality  follows  from  <■*).  Our  estimate  comes  from  this 
last  expression.  Observe  that  for  fixed  v,  the  maximum  of  tarti'/2  will 
occur  when  -©=3->f/?,  i.e.,  where  the  situation  is  as  in  figure  5. 


Figure  5. 


Thus,  jtanSf/2j  <v/2y,  so  that 
3 

lEl*T^T  <2si"l*™|-  7T  >- 

+  /  v 

Also,  for  small  enough  ■©,  me  certainly  have  |sin3j|<v/y,  so  the  final 
estimate  for  this  Fresnel  Approximation  error  becomes 
|e|  i  3v-V<8y=»3. 

We  now  obtain  a  similar  estimate  for  the  near-normal  Fresnel 
Approximation.  Again  referring  to  figure  3  above,  this  approximation 
is 

R-rX<  <x+v)a-xa)/2ynsxv/y  +  va/2y. 

It  arises  by  differencing  the  truncated  power  series  expansions  of  the 
square  roots  in 

R=y Jl+< <x-*-v)/y)a,  r=y Jl+( x/y )*. 

Since  near-normal  incidence  is  being  considered  here,  to  estimate  the 
error  we  begin  by  assuming  that  [xj<v.  An  exact  equation  for  R-r, 
also  given  above,  is  (assuming  -©>0 ) 

R-r=xv/r  +  ( va-( R-r >a) /2r . 

Thus,  the  error  which  we  wish  to  estimate,  E*,  is 

E*=xv< 1/y-l/r )+va< 1/y-l/r >/2+<R-r )a/2r 
Note  that  E*>0,  so  that  the  exact  value  of  R-r  is  less  than  the 


approximate  value. 


R-r<xv/y+v2/2y<3v2/2y .  Also, 
r-y=x2/ ( r+y ) <xa/2y . 


Putting  these  estimates  into  the  above  expression  for  E* 

yields 

0<E*<x3v/2rys  +  x2v2/4ry2  +  9v"V8ry2. 

Since  x<v  and  r>y, 

0<E*<15v*/8u*. 
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The  presence  of  errors  resulting  from  such  causes  as  clutter 
scatterer  returns,  element  thermal  noise,  electronic  counter'" measures, 
antenna  location  errors,  and  multipath  effects  often  gives  rise  to 
instability  in  nonlinear  high-resolution  beamforming  techniques.  One 
reason  is  that  these  methods  can  place  undue  emphasis  upon  unstable 
eigenvectors  of  the  noise  covariance  matrix  Clll.  Hence  the  question 
of  how  errors  in  the  entries  of  the  noise  covariance  matrix  effect  the 
basic  properties  of  an  array,  such  as  its  gain,  is  of  interest. 
Preliminary  results  in  this  regard  may  be  found  in  Byrnes  and  Sullivan 
C193,  and  a  continuation  and  expansion  of  these  ideas  is  presented  in 
Appendix  F. 

Beamforming  enables  an  array  to  act  as  a  spatial  filter,  by 
enhancing  detectability,  resolution,  and  directional  measurement  of 
plane  wave  signals.  There  are  certain  applications,  however,  where 
beamforming  can  offer  only  peripheral  assistance,  and  other  methods 
are  required.  An  important  example  in  this  regard  is  low-frequency 
detection  in  the  ocean.  Propagation  loss  modeling  is  a  fundamental 
tool  applied  to  the  detection  of  targets  at  medium  to  long  ranges,  and 
the  Parabolic  Equation  Method  (PEM)  has  recently  become  a  basic 
technique  in  the  attempt  to  improve  the  solvability  and  accuracy  of 
such  models.  Our  investigations  into  the  PEM,  including  the  discovery 
of  certain  inaccuracies  in  its  derivation,  are  described  in  Appendix 
C. 
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VI.  ROBUSTNESS 


We  now  report  on  additional  progress  in  determining  the 
robustness  of  trigonometr ic  polynomials  with  unimodular  coefficients. 
As  indicated  in  C18  3,  we  do  this  by  estimating  the  maximum  of 
tr igonometric  polynomials  with  random  coefficients.  These  random 
polynomials  represent  the  deviation  in  frequency  response  of  a  linear, 
equispaced  antenna  array  caused  by  coefficient  inaccuracy  in  the  basic 
unimodular  polynomial  of  the  array. 

Let  us  consider  this  exposition  the  random  cosine 

polynomial 

P<t,w)  =  X  *  ( w )  cos  k t , 

Lil 

where  the  Xt«are  independent,  standard  normal  random  variables,  denoted 
Xu  i. i.d  N(0,1>.  Our  results  apply  immediately  to  the  sine  polynomial 

XfccSin  kt  and  may  also  be  generalized  to  the  random  tr igonometric 

x-\ 

yN 

polynomial  £XueA**. 

-V> 

We  define 

M(w>  =  sup  |P(t,w)j 

and  we  wish  to  estimate  ,  <M<  c<a)’.  Nate  that  n,  the  degree  of 

P(t,w),  is  a  fixed,  positive  integer.  Simple  upper  bounds  for  P'{M>  <x}- 
were  obtained  in  our  earlier  work,  and  we  now  seek  a  lower  bound  for 
this  probability  to  complete  our  estimation.  In  the  future,  we  hope  to 
expand  this  work  to  include  more  general  random  coefficients.  For  now, 
we  limit  ourselves  to  the  normal  case. 


We  give  two  results.  The  first  is  less  precise  and  loss  useful 
than  the  second,  but  we  include  it  to  introduce  a  method  of  promising 
simplicity  which  we  hope  to  exploit  later.  The  second  is  a 


case. 


generalization  of  a  result  of  Salem  and  Zygmund  C 38 1  to  our 
After  stating  and  proving  these  results,  tue  shall  compare  their 
precision  with  other  similar  estimates. 

a  «) 

p-  _ _ 

Proposition  Is  (1)  PCM>j«>  >  1  -  1 _ where  0<>S</n/2, 

or  equivalently 

ci-£) 

(2)  PCMC/&X  ~r=  where  0</5  <  \fn72. 

~  vn/x 


Proposition  2s  Let  &€ < 0, 1 ) .  Then 


- - io.  &  ' 

( 1 )  PCM>  1/n  logn  +  -A- 


In  particular,  if  6  =  ( logn )_1 , n^3,  we  have 

-»  * 

C2)  PC.H>  Vn  logn  -  ^  ?  ‘‘  ^ 

—  *~f  'ey,  ^  —  t  4-  - -  vA 

2.^ 


Proof  of  Proposition  Is 

We  prove  Proposition  1  (2), 


i  ^  a 

Let  >  >0.  PCM</3  >=PC-y  M=*>  /3  >e  ^  >. 

By  Markov's  Inequality, 

*>»  7) 

__  m 

PCM</3  l<e  ^  E(e  x  ). 

Evidently,  M*<w)  >Pa(  t, w)  for  every  t<£0,2TTl,  and  since 

A 

P(t,w)~N<0,  cos*kt),  we  have 
K-l 


-45- 


^  v'  N*  v  '  *v  v  v  v  V  V>  v  V  ->  *>  \»  •«>.%»  \>  -V.  «>  •%>  \N  \\\y  -VI  - 


p  ( 


( ^  X2  (chi-square  with  or.e  degree  of  freedom). 

;  l  ^ 

Therefore,  using  the  fact  that  the  moment  generating  function  <J)  of  a 
X2  random  variable  is  <}>(■£  )  =  \/^\ '■>i>  for  \<i/2. 

-i*  -2s p*  /.  -3  rr  xl  x 

E<e  -m)^E<e*p  >=E(e  *•  u-*  '  ) 


( 10) 


\rnrgf  «,.v 


Combining  (9)  and  (1(8)  yields 


PCMOS  X 


e 


\ji 


Now  let  us  find  >B  and  t<C0,21T3  which  minimize  the  right  side  of 
the  above  inequality.  Clearly,  t=0  is  the  correct  value  giving 

2_cosak*0=n  and  thus 

Vn\ 


A  & 


PCM<>S  K 


=»h(7s). 


v-) 


(11) 


To  find  7*  which  minimizes  h^’  ),  we  have 


h‘  <7  )= 


|/3  r  <L  O  +  ~  -O'  -  \ 
*  L-a  x  ,i  3 


(,»  -V-  —  — 


, v"  J 
o. 


j  V  ^ 

Thus,  h  >=0  <==£>  •—  ( 1+T  n)~  --  =0  «=s>  »-4 

*  1  S  ^ 


and  ?*  is  obviously  a  minimum  point  of  h.  To  insure  ^  >0,  we  require  ft 

\j~^  and  substituting  the  critical  value  of  ?)  in  equation  (11), 

*»• 

have 


we 


\T"  / 


giving  the  result. 
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Proof  of  Proposition  2: 


We  follow  the  method  in  Salem  and  Zygmund  C38J,  modifying 
their  results  to  fit  our  case.  We  begin  mi  h  a  preliminary  lemma. 

Lemma  Is  Let  Q(w)  be  a  non-negative  random  variable  satisfying 
E(Q)>A>0  and  E(Qa)<B  where  E  denotes  the  mathematical  expectation. 
Then,  if  cfe(B,  1>,  P<Q>ck)> ( l-<£) aA2/B. 
proof 

Put  JB&{w:Q(w) >dA^ .  Then  C  Q(w)dP(w)<c^AP(jyt!;)  <t^A 

#>t 


^  QdP  ■  ^„QdP  -  jQ  QdP  >  EQ-^A  >A-t$A  =  (l-cf)A. 


Also,  ^QdP  ^  ladP)  *"  (^QadP) 

<  P(^T)  '/agi/a. 


u 


Combining  the  above  inequalities,  w©  have  (  l-<£)  A<P(^)  *  which 

yields  the  conclusion  of  the  lemma. 

We  wish  to  apply  lemma  1  to  the  random  variable 


o 


To  do  this,  we  must  verify  that  ? to)  satisfies  the  hypotheses  of  the 
lemma,  and  we  must  obtain  the  constants  A  and  B. 

Bm  Tonelli’s  theorem^ 


E<  I ) 


--f 


*  PC* 


E(e 


)dt  «— > 


nr  p, 


e  °  * 


dt 
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m 


ft 

4 


f 


using  the  form  of  the  moment  generating  function  of  P(t,w)  already 


obtained.  With  the  simple  identity  cos2kt  =— C l+cos2kt ],  me  find 


E(I)  = 


■  n 


2-  O 

2L  ^cOS9-W- 


*4  v_c_  \ 


air 


>  s'1 


-XT 


ATT 


^  ( 1+  "3~2Lc:os2kt  )dt  =  e  g 


So  we  have  Ed)  >  e  n  =  A  in  lemma  1. 

We  now  require  a  number  B  such  that  E(Ia)<B  to  satisfy  the 
second  hypothesis  of  lemma  1.  To  do  this,  we  write 


\  ^  jcr 

_  /  /*  /*l4- 


I3(w)  =— -5r  j"  y  e  "  U  ^ 


1  /  /"Ujrv.  w  ,v>u)j\ 

Then  Ed3)  *  j  E(e  jdt  dfc 

o  o 

using  Tonelli’s  theorem. 


P(t,w)  +  P(T,  w)  =  £L  X*(w)  (cos  kt  +  cos  k“f ) 


V.--I 


^  N(0,  Z.  (cos  kt  +  cos  k"id3) 

*-) 

because  the  sum  of  independent  normal  randan  variables  is  normal.  We 


thus  have  that  the  moment  generating  function  of  P(t,w)  +  Pft^w)  is 


(j 


Ele_ 


-  e.  't 


Then 


Jir  *tr  ^  £l  («sbVt*<o',V-r} 

^  f  (  9-  >t-' 

Ed3)  » — *  lie  dt  dt: 

MTT  o  *o 

Now  (cos  kt  +  cos  tcf)  a=  l+^cos2kt  cos2k't"  +2cos  kt  cos  Wt. 
Put  S  =  S(t,'t)  =  Zi~  cos2kt  +  ~cos2kf  +  2cos  kt  cos  ktj • 
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f- 


This  gives 


E(I2)  = 


Asr  ?T 

e  •*- 


dt  df. 


To  bound  the  final  integral,  we  mill  use  lemma  4.2.2  of  Salem  and 


Zygmund  which  we  state  for  functions  of  two  variables. 

Lemma  2  < Salem  and  Zygmund) 

Let  g<x,y),  a<x<b,  c<y<d,  be  a  bounded  real  function.  Suppose 
that  |g(x,y)|<A  and  ( b-a)-1 < d— c) -t  ^  y  g2Cx,y)dxdy  *  B.  Then,  for  any 
positive  number  u, 

d  V> 

<  b-a )-1  (  d-c )  “  1  ^  ^  e««<M.w>dxdy  <  1 +u/B+ Jj; e“A 


In  the  case  when 


replaced  by 


V  « 

1  (. 


|(x,y)dxdy  ~  0,  this  inequality  can  be 


(b-a)-Md 


!— c)~*  ^  {  e*^0<M» 


M,dxdy  <,  1+^*-  e“A. 


We  apply  lemma  2  to  the  function  S=S(t,'Y)  for  0<t<2lT, 
IT.  To  verify  the  necessary  hypotheses,  first  note  that 


>TT  atr 

5  {  Sd^JdtdiT  *  0. 


The  set  of  3n  functions  Ccos2kt,  cos2k1T,  cos2kt  cos2kt  :k=l,*n]  is  an 
orthogonal  system  over  the  square  of  integration.  This  leads  to 


^  $  Sa(t,T)dt  c ft 
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air  q/tr 


=  tirp"  C  (  (z  cos22kt+-q  cos22k;£+4cos2kt  cos^kt^)  dtdy 
V.r|  /o 


=<^  +j  +l)n  =  ^  n  =  B  of  lemma  2. 


Evidently,  [S(t,‘£)|  <  3n=A  of  lemma  2,  Thus,  the  hypotheses  of  lemma  2 
are  satisfied  and  ute  obtain  ,,  with  u=  , 


Applying  this  to  equation  (12)  gives 


E<I2)  <  e ^  C  1+rr;  e*  " 3  =  B  of  lemma  1, 


Recall  that  we  have  obtained  E(I)  >  s'4  =  A  of  lemma  1.  Thus,  we 

apply  lemma  1  to  obtain 


<. '  —  <5  ) 


PO  >  >  >  7T~r 


for  c£<(0,  1  ) . 


To  prove  proposition  2,  ute  note  that  for  every  /  >0,  t€X’0,2Ti3 
and  w,  we  have  e  <  e  Applying  )  dt  to  both  sides  of  the 
inequality  yields  I(w)<e  for  each  r  >0.  If  we  use  this  with  (13), 


we  get 


log<i  ^  n)-  =  p(e  >  6  eM  ). 


(.  \ -6  j 


>  P{l  >  6  e  )  >  — - 5 - 5: 

1  +—*/«_  A 
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Now  put  T'  =4  y^in  the  inequality  to  obtain  proposition  2(1)  and 
then  substitute  <£  =<logn)_1,  n>3,  to  obtain  proposition  2(2). 


Discussion  of  Results 

We  now  examine  propositions  i  and  2  and  compare  them  to  other 
similar  estimates.  An  excellent  standard  of  comparison  is  the  recent 
asymptotic  result  of  Turkman  and  Walker  C443  that 


tor  every  real  number  x.  Roughly,  this  says  that  for  large  n,  M=M„  is 
centered  near  s/rTTogn  with  dispersion  of  the  order  of  v .  The 
difficulty  in  applying  this  result  is  that,  for  a  given  n,  it  is  not 
known  how  close  the  distribution  of 


lo« 


■<V\  ~  Vr>  n  ~ 


r» 


t  / 


'<><<  V\ 


is  to  the  limiting  distribution. 

Comparing  this  result  to  Proposition  2(2)  shows  that  our 
centering  constants  closely  resemble  those  of  Turkman  and  Walker  and 
suggests  that  the  bounds  obtained  in  Proposition  2  are  indeed  tight. 
In  addition,  the  lower  bound  of  Proposition  2(2)  is  better  than  a 
similar  bound  obtained  by  Gersho,  Gopinath  and  Odlyzko  C221  in  the 
second  part  of  their  theorem.  The  Gersho  bound,  however,  applies  to  a 
more  general  class  of  random  tr igcnometr ic  polynomials.  It  is  obvious 
that  the  bounds  of  Proposition  2  are  much  superior  to  those  of 
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Proposition  1,  which  applies  only  when  JB<  i/n/2,  which  is  far  below  the 
central  values  of  M,.  Proposition  1  has  been  included  only  to 
illustrate  a  simple  but  promising  method  which  we  hope  to  exploit  in 
the  future. 

The  future  program  is  to  generalize  the  methods  of  Proposition 
2  to  trigonometric  polynomials  with  independent  random  coefficients 
which  are  not  necessarily  normally  distributed.  Of  particular  interest 
is  the  case  of  uniformly  distributed  coefficients  which  represent 
digital  or  round-off  errors.  While  the  applicability  of  the  methods  of 
Salem  and  Zygmund  to  the  attainment  of  an  upper  bound  for  P(M  > c<)  has 
been  noted  by  Gersho,  et  al  an'd  others,  nobne,  to  our  knowledge,  has 
used  these  methods  to  obtain  a  lower  bound  as  we  have  done  in 
Proposition  2.  This  development  offers  hope  for  improved  probablistic 
estimates  of  the  maximum  of  trigonometric  polynomials  with  general 
random  coefficients. 
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Tuscany,  Italy,  18  August,  1987 

5.  University  of  Montreal  Mathematics  Research  Institute,  10 
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I.  Invited  Paper  Presented  by  Stephen  Boyd 

1.  Naval  Weapons  Center,  China  Lakes,  CA,  20  April,  1987 

II.  Consultations  on  Potential  Air  Force  and  Navy  Applications  of 

Prometheus  Ideas. 


1.  Stephen  Boyd  and  James  Byrnes  consulted  with  Bob  Dinger, 
Bill  Altop,  Gary  Heuier,  Francis  Canning,  and  Mike  Mumford  at 
the  Naval  Weapons  Center,  China  Lakes,  CA,  20  April,  1987. 

2.  James  Byrnes  consulted  with  Emanuel  Vegh  and  Bill  Gabriel 
at  the  Naval  Research  Laboratory,  Washington,  DC,  24  April, 
1987. 


3.  James  Byrnes  consulted  with  Neil  Gerr  and  Jim  Smith  at  the 
Office  of  Naval  Research,  Arlington,  VA,  24  April,  1987. 

4.  James  Byrnes  consulted  with  George  Behnke  of  the  NITRE 
Corporation,  27  April,  1987  (this  meeting  took  place  at  the 
University  of  Maryland). 

5.  James  Byrnes  consulted  with  H.  S.  Shapiro  at  the  Royal 
Institute  of  Technology,  Stockholm,  24-26  June,  1987. 

6.  James  Byrnes  consulted  with  numerous  people,  especially  E. 
K.  Miller  of  the  Rockwell  Science  Center  and  the  University  of 
Kansas,  N.  H.  Farhat  of  the  University  of  Pennsylvania,  Y, 
Rahmat-Samii  of  the  Jet  Propulsion  Laboratory,  Tapan  Sarkar  of 
Syracuse  University,  and  Anton  G.  Tijhuis  of  Technische 
Hogeschool  Delft  at  the  NATO  Advanced  Study  Institute  on 
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Electromagnetic  Modeling  and  Measurements  for  Analysis  and  for 
Synthesis  Problems,  Tuscany,  Italy,  10-21  August,  1987. 

7.  James  Byrnes  and  Donald  Neu/man  consulted  with  Bob 
Mailloux,  Hans  Steyskal  and  Jeff  Herd  at  the  Rome  Air 
Development  Center,  Heinscomb  AFB,  Bedford,  MA,  24  August, 
1987. 


8.  James  Byrnes  consulted  with  Jack  Wallace,  Mark  Godino  and 
Mike  Carpenter  at  Wright  Patterson  AFB,  Dayton,  OH,  25  August, 
1987.  This  meeting  was  also  attended  by  Arje  Nachman  of 
AFOSR. 

9.  James  Byrnes  consulted  with  Stanley  Chamberlain  and  E.  Ted 
Bick  at  Raytheon  Submarine  Signal  Division,  Portsmouth,  RI,  22 
September,  1987. 

10.  James  Byrnes  consulted  with  many  SBIR  officials,  DOD 
representati ves,  and  representati ves  of  several  large  defense 
contractors  at  the  High  Technology  Conference,  Salt  Lake  City, 
UT,  17-18  November,  1987. 

1  j.  .  James  Byrnes  visited  Sun  Microsystems  in  Mountain  View, 
CA,  to  evaluate  the  appl icabi 1 i ty  of  their  engineering  work 
stations  to  the  Prometheus  Inc,  Phase  II  effort,  to 
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The  L*  Norm  of  a  Polynomial  With  Coefficients  +1 

Donald  J.  Neuman  and  J.  S.  Byrnes 
Prometheus  Inc. 

103  Mansfield  Street 
Sharon,  MA  020o7 

A  classic  unresolved  question  regarding  polynomials 

P(z)=  <*<n)z“ 

k'-O 

with  coefficients  <  „  >  =<£*= +1  is  whether  the  maximum  modulus  of  such  a 
polynomial  on  the  unit  circle  can  be  nx''a+o(nly'a) .  As  shown  by  Kahane 
C3  3,  if  complex  coefficients  of  modulus  1  are  allowed  then  not  only  is 
it  possible  for  this  property  to  be  satisfied,  but  the  minimum  modulus 
can  be  n1''a+o(nly'a)  as  well.  Specifically,  Kahane  proved  that  for  any  n 
there  is  a  polynomial  of  degree  n  with  coefficients  of  modulus  one  whose 
modulus  everywhere  on  the  unit  circle  is  n1''a+0(n3<'**log  n). 

Erdos  C 2 3  had  conjectured  the  existence  of  a  c>0  such  that,  for 
any  polynomial  P  of  the  types  described,  ||p|j  >  ( l+c)nW5*„  Clearly  the 

cO 

Kahane  result  disproved  this  conjecture  for  the  modulus  1  case,  but  the 
situation  for  coefficients  +1  remains  open.  Employing  an  elegant 
construction  Shapiro  C6,53  demonstrated  the  ach ievabi 1 i ty  of  the  order 
of  magnitude  n*'"2,  but  the  maximum  modulus  of  the  Shapiro  polynomials  is 
about  <2n)1''2. 

Since  the  L2  norm  of  any  such  polynomial  is  n1'2,  by  the 

Parseval  Theorem,  and  the  Lp  norm  for  any  p<C*3  is  a  lower  bound  for  the 

maximum  modulus, it  is  natural  to  look  at  the  Lp  norm  for  some  p>2. 

Motivated  by  these  considerations,  as  well  as  the  inherent  traetability 

of  the  L*  norm,  we  examine  the  L*  norm  of  such  polynomials.  As  one 

might  expect,  this  leads  to  several  interesting  combinatorial  questions. 

We  provide  answers  to  some  of  these,  and  conclude  with  a  refined  version 
of  the  Erdos  conjecture. 


Throughout  the  paper  n  mill  be  a  positive  integer,  P(z)  will 
denote  the  previously  indicated  polynomial  of  degree  n-1  with 
coefficients  +1,  and  z  will  lie  on  the  unit  circle,  0<<X1. 
All  integrals  will  be  over  -$<110, 13.  We  begin  with  a  Lemma. 


Lemma 


Prgof 


L* 


r  C-IL  L  . 


=■  ^constant  terms. 

Since  a  constant  term  occurs  in  this  product  if  and  only  if  j+k»l+ro, 
the  result  follows  immediately. 

.  ^ 

Of  interest  is  the  expected  value  E(  (fp([  ),  if  the 
coefficients  £j  are  chosen  at  random. 

Ibesnsa  i  etflp!L)«2na-n. 

u 

PLSaf  Clearly  if  exactly  3  of  the  indices  j,k,l,ra  are 
identical,  or  if  at  least  3  of  them  are  different,  then  E(£.*£ 

r  )  ^ 

'V520*  lt  Prefer©  follows  from  the  lemma  that 


u> 


Efiipf 

'  ^  ^ 

0t  j-M  IrJ  JU£ 

For  each  of  the  n(n-l)/2  pairs  of  integers  p,q,  0^p<q<n, 


L  .<<L  C  ' 
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there  are  4  terms  appearing  in  ( 1  > ,  being  typical.  The  only 

other  terms  in  <1)  are  ,  0<p<n.  Since  all  of  these  terms  equal  1, 

¥-  ' 

E<  IIpII  ^)=n+4(n>  (n-1  )/2=2n*-n, 

L 

completing  the  proof  of  Theorem  1. 

We  now  observe  the  improvement  that  is  achieved  when  this  random 
choice  is  replaced  by  the  Shapiro  coefficients.  Shapiro’s  polynomials 
are  defined,  together  with  his  auxiliary  polynomials  Q,  by  the 
recurrence  formulas 


P* ( z ) =Qm ( z ) - 1 , 

<2>=P«,(z)+z2rTQm(z) ,  (2) 

n 

Q«,^.».(i:)=Pm(2)-2a  Q„(z),  m>0. 

For  the  interested  reader  we  point  out  articles  by  Brillhart  et  al  C 1 3 
and  Mendes-France  and  Tenenbaum  C43.  The  former  shows  that  the  Shapiro 
coefficients  can  be  defined  directly  from  the  binary  representat ion  of 
the  order  of  the  polynomial,  while  the  latter  relates  them  to  paper- 
folding  sequences. 


Theorem  2  If  n=2K  and  P(z)  is  the  Shapiro  polynomial  of  degree 


n— 1,  then 


||p||^.=  (4na-<-l)**n>/3. 

Ul 


Proof  It  follows  directly  from  <2>  that 

|p— *<z)la  +|o«**(z)f  a=2(  lp«(z)|  *  +  l0«.(z)J  *) 


so  that,  as  observed  by  Shapiro, 


|P«(z)|  *-HQ«.<z)J  *s2r~*. 


Now  (2)  and  <3)  uield 


jp*.*t<z>|  *=2~* *  +  2Re( z*  Qm(z>Pmlz)>, 
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VK**  p"  is  composed  solely  of 

Next  we  observe  that  *  M  «n  »i 

.  __  +hal:  it  can  be  thought 

frequencies  which  are  positive  powers  of  2,  so 

of  as  Q  ? .  -her.  F  is  the  -reversed-  polynomial  of  P  -  Thus 

siuK-*mn  *  *  ■ 

Since  fsoF  is  analytic  and  0  at  the  origin. 


so  that 


I  »  —  —  - - 

0--  /f*  a  &(frX)  «  f(&. if-  Jk*)\ 
"  f(M)1  -- 1 

_  K 


Altogether  then,  we  have 

I*  .  .  m  f  IP  I1/. 


ivj  -t“*-  tfifj'i-r-itj) 

-  jl'h+3  J Ifinl* 

The  remainder  of  the  "roof  i.  now  simply  on  k.  The 

„  _  ,  _  p  (.)ai,  Furthermore,  from 

result  is  obviously  true  for  k-B,  •*«»  P* 

(5)  and  the  inductive  hypothesis 

■  .  L 


T.k+1 


-  (-1) 


xfcfl  1 


:lt+2  -  f-i)* 


it  follows  that 

rn  it  lfcrt  .zi"  ■  J 

fa, I  »  2  -2 - 3 -  = - 3 - -  •  "  T"*f- 

This  completes  the  proof  of  Theorem  2. 

o  Unlies  that  the  L*  norm  of  the  n-i“* 
Note  that  Theorem  2  implies 

degree  Shapiro  polynomial  is  asymptotic  to  *  Uses  the  fourth  root 
Of  */3-l.  07437^.  Based  upon  extensive  numerical  evidence  eopioyin, 
the  Bose-Einstein  statistic,  methodology  of  St.tisticai  Hechanics.  w. 
conjecture  that  the  Shapiro  polynomial.  do  not  give  the  «.i—  ^ 
„orm  among  aii  polynomial.  of  the  same  degree  with  coefficient.  ±», 
but  that  this  minimum  L»  norm  is  asymptotically  Vn  times  the  fou 
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root  of  6/5**l .  04664)/"n.  Observe  that  the  truth  of  this  conjecture  would 
imply  that  of  the  Erdos  conjecture  mentioned  earlier,  with  c» < 6/5 ) * 
l  a*.  04664. 
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An  Analysis  of  the  Parabolic  Equation  Method 
and  its  Applications 
J.  S.  Byrnes 
Prometheus  Inc. 

Abstract 

The  parabolic  equation  (PE)  method  of  solving  the  elliptic 
wave  equation  in  underwater  acoustics  was  developed  to  handle  the  case 
where  the  index  of  refraction  n  is  dependent  upon  range  as  well  as 
upon  depth,  and  possibly  azimuth.  It  is  used  extensively  in  low- 
frequency  propagation  loss  modeling.  Fundamental  to  the  derivation  of 
all  acoustic  PE’s  is  the  assumption  that  a  certain  commutator  is 
negligibly  small,  and  may  be  ignored.  One  purpose  of  this  paper  is  to 
point  out  that  the  vanishing  of  the  commutator  actually  implies  that  n 
is  range- indepen dent.  Thus,  in  many  cases,  the  original  reason  for 
the  PE  method  is  defeated  when  this  commutator  is  assumed  to  be 
negligibly  small.  The  second  purpose  of  this  paper  is  to  describe  and 
analyze  the  current  state  of  the  art  in  the  low-frequency  applications 
of  the  PE  method,  and  to  compare  it  to  other  methods  of  solving  the 
acoustic  wave  equation. 

Key  words:  Parabolic  Equation  Method,  Propagation  Loss  Modeling 
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An  Analysis  of  the  Parabolic  Equation  Method 
and  its  Applications 
J.  B.  Byrnes 

Prometheus  Inc.,  103  Mansfield  St.,  Sharon,  MA  02067  and 
University  of  Massachusetts  at  Boston,  Boston,  MA  02125 

I.  Introduction 

Knowledge  of  how  sound  propagates  in  the  ocean,  for  different 
environmental  conditions  and  source/receiver  configurations,  is 
necessary  in  order  to  achieve  optimum  sonar  design.  The  estimation  of 
the  spatial  properties  of  the  sound  pressure  field,  as  a  function  of 
source  frequency,  is  the  ga,*l  of  ocean  acoustic  modeling.  To  perform 
this  modeling  effectively,  the  acoustic  loss  mechanisms  in  the  ocean 
must  be  accurately  handled.  These  loss  mechanisms  include,  in 
addition  to  geometrical  spreading  loss,  bottom  reflections,  volume 
absorption,  and  scattering. 

Various  methods  of  handling  acoustic  wave  propagation  problems 
exist,  because  no  individual  method  adequately  deals  with  the  many 
types  of  problems  and  conditions  that  one  is  likely  to  encounter  in  a 
variable  ocean  environment.  These  methods  include  ray  theory, 
fastfield  theory,  normal  mode  analysis,  and  tne  parabolic  equation 
method.  An  important  test  of  any  particular  method  is  its  ability  to 
handle  substantial  variations  in  environmental  parameters,  and  to 
accomplish  this  while  using  only  a  reasonable  amount  of  computer  time. 

The  acoustic  wav©  equation  cannot  be  separated  if  the 
environment  varies  in  both  depth  and  range,  so  that  direct  numerical 
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integration  is  required.  Since  the  basic  equation  represents  a  three- 
dimensional  boundary  value  problem,  there  are  currently  no  practical 
methods  available  to  solve  the  problem  in  this  way.  Directly 
integrating  the  acoustic  uave  equation,  with  appropriate  boundary 
conditions,  remains  impractical  due  to  the  excessive  computation 
times,  even  taking  into  account  the  recent  leaps  forward  in  computer 
technology.  Thus  some  simplifying  assumptions  are  always  introduced 
before  the  wave  equation  is  solved.  The  parabolic  equation  (PE) 
method  is  thought  to  represent  a  viable  alternative,  as  it  is  an 
approximation  to  the  uave  equation  which  lends  itself  to  practical 
numerical  solution. 

II.  Range  Independence  in  the  Parabolic  Equation  Method 

Study  of  the  parabolic  equation  (PE)  method  was  initiated  by 
Tappert  and  Hardin  C42  3,  in  an  attempt  to  deal  with  an  ocean 
environment  in  which  the  sound  speed,  and  hence  the  index  of 
refraction,  is  dependent  upon  range  as  well  as  upon  depth,  and 
possibly  azimuth.  Since  1974  the  PE  method  has  become  the  standard 
way  of  dealing  with  this  ’range-dependent"  case.  However,  this  is 
theoretically  incorrect,  as  a  basic  assumption  necessary  to  the 
derivation  of  the  PE  actually  implies  that  the  index  of  refraction  is 
independent  Of  range.  To  see  why  this  is  so,  it  is  necessary  to 
examine  the  derivation  of  the  PE  in  detail.  A  convenient  reference  in 


this  regard  is  £38 3,  and  we  begin  with  equation  (8)  of  this  paper. 

+2iK»  5^“  ♦k**(n*-l>u«0.  £  38 . 8  3 

Here,  r,  z,  and  -O' are  the  range,  depth,  and  azimuthal  variables 


respectively,  u=u(r,z,0’)  is  the  (unknown)  modulation  of  the  spatially 


varying  part  of  the  acoustic  pressure,  c=c(r,z,^)  is  the  actual  sound 
speed,  is  a  reference  sound  speed,  n=n  (r ,  z,ff>=Cg /c  is  the  index  of 
refraction,  w  is  the  acoustic  frequency  in  rad/sec,  and  ka=w/c*  is  the 
reference  wave  number.  As  usual,,  the  farfield  assumption,  k«,r  >>  1,  is 
made. 


Equation  C38.B3  may  also  be  expressed  in  operator  format: 

(  fV+2ik«f-  +  -fr  +4  -sL?  +ka  Mnv- 1 )  )u=0.  C38.ll  3 

d<  J*  <■ 

Two  seemingly  different  methods  of  deriving  the  fundamental  PE 
have  appeared  in  the  literature,  one  employing  C38.83  and  the  other 
C38.il  3.  To  derive  the  PE  from  C38.83,  the  "paraxial  approximation," 
u, rXfl,  is  made.  One  supposed  justification  for  this  is  the  assumption 
that  u  varies  "slowly  as  a  function  of  r  on  a  wavelength  scale,  i.e. , 
JUr}<<3K,uJ,  which  implies  that  |urrj<<2ik<,u»j. "  C29,  p.475  3.  This 

implication  is  clearly  incorrect,  yet  it  is  the  foundation  upon  which 
the  paraxial  approximation  is  based. 

The  second  standard  method  of  obtaining  the  PE  is  to  replace 
C38. 11 3  with: 

( +  ika  -ik<,Q)  (  -*‘ik(9-* ika Q)u=<3,  E 38*  12  3 

Ov  <K 

where  the  operator  Q  is  defined  by 


Q*(l+!n‘-l)+pfr  ♦rrr  >*.  C 38. 13  3 

The  difference  between  the  operators  in  C 38. 113  and  C 38. 123  is 
the  commutator  C^D^Q—QD^,  where  Dr  indicates  the  operator  of  partial 
di f ferentiation  with  respect  to  r.  Thus,  in  order  to  obtain  C 38. 123 
from  £38.113,  it  is  necessary  to  assume  that  C  vanishes  (i.e.,  is 
negligibly  small).  After  this  is  done,  the  second  factor  in  C 38. 123, 
which  represents  the  incoming  wave,  is  neglected,  yielding 
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C  38.  14  3 


Jo- 

— +ik*u=ikaQu. 

When  the  square  root  operator  Q  in  C 38. 14 3  is  replaced  by  the 
linear  terms  of  its  Taylor  series,  the  resulting  equation  is  identical 
to  that  obtained  by  applying  the  paraxial  approximation  to  C38.83, 
namely  the  standard  three-dimensional  (3D)  PE, 


if  . . . 


A. 


C  38 .93 


To  see  that  the  two  basic  assumptions  employed  as  alternative 
means  of  deriving  the  PEf  namely  the  paraxial  approximation  and  the 
vanishing  of  the  commutator  C,  are  in  fact  the  same,  assuming  that  the 
incoming  utave  is  ignored,  begin  with  the  equivalent  equations  C38.83  and 
C  38. 113.  Then  s 


|  ,  .  i _  j 

<ss>(^7'  +  ik<,-ik<»Q)u?*8<=s>(  ^  +ikd-ik„ Q)  +ika-fikd<3>ut»0 

<=>Cup43. 

Therefore,  it  is  necessary  to  assume  that  Cf=0  in  order  to  derive 
the  PE.  However,  in  the  azimuthally  independent  two-dimensional  (2D) 
case  the  vanishing  of  the  commutator  actually  implies  that  n  is  range- 
independent,  while  in  the  3D  case  the  assumption  means  that  n  is  ’very 
nearly*  range- independent.  Thus,  in  many  cases,  the  basic  theoretical 
reason  for  the  PE  method  is  defeated  when  the  commutator  is  assumed  to 
be  negligibly  small. 
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The  computation  required  to  understand  this  is  straightf orujard. 


First,  observe  that  the  commutivity  of  D  with  Q  implies  that  of  D  with 
Q*,  since  DQa=  DQQ=QDQ=QQD=QaD .  Hence,  the  (approximate)  vanishing  of 
DvQ-QD^  implies  that  of  D^Q^Q*  D,..  However, 

CDvQa  -  QaD„ 3  u  =  C ( 1+X+Y)  -  (1+X+Y)D,.3  u 
=  2 ( n n ,»  —  kD~2  r-3  a )  u ,  where 

X=na-l+k&-aD,>,  Y=(kar)-aD*a. 

In  the  2D  case  D^=0,  so  Cu=0  implies  that  nv=0,  i.e.,  the  refractive 

index  is  range-independent.  In  the  3D  case,  Cu-0  means  that 
nv.=!Da.a(u)/nk3  *-  r1.  Combined  with  the  farfield  assumption,  it  is  clear 
that  n,.*®„ 

It  is  appropriate  to  observe  that,  in  their  discussion  of  the 
approximation  resulting  from  neglecting  the  commutator,  Seigmann,  et  al 
might  have  reached  the  same  conclusion.  A  basic  consideration  at  this 
point  in  C383  is  that  the  term  k0-1  (n2~l)-‘n^  must  ba  negligible.  Since 
ni=4  (and  in  fact  n*l  at  various  points),  this  requires  that  n^p0.  The 
primary  error  occurs  in  equation  (34)  of  C38  3,  where  “nan*  Should  be 
replaced  by  “min,"  and  min  J —  l J  e0. 

Thus,  it  is  necessary  to  consider  the  PE  method,  from  a 
theoretical  viewpoint,  as  an  alternative  way  of  solving  the  range- 
independent  case,  where  it  is  known  to  perform  efficiently  (see  C  6  3,  for 
example),  as  opposed  to  a  method  of  solving  the  range-dependent  case. 
For  the  purpose  of  practical  numerical  computation,  the  range-dependent 
case  has  been  attacked  by  approximating  n  with  a  function  that  is 
piecewise  constant  in  r.  Numerical  schemes  employing  the  PE  method, 
such  as  the  wide-angle  split-step  Fourier  method  C433  and  the  wide-angle 
implicit  finite  difference  (IFD)  scheme  C53,  incorporate 


this  idea. 


III.  The  Standard  PE  Method,  and  Alternatives 

There  are  several  basic  assumptions  which  must  be  made  in 
order  to  apply  the  standard  PE  method.  These  includes  the  fai — field 
approximation;  the  index  of  refraction  is  piecewise  constant;  there  is 
only  an  outgoing  wave;  there  is  no  scattering;  the  frequency  is  less 
than  500  Hz;  the  propagation  angle  is  small  (less  than  about  15*). 
One  way  of  possibly  avoiding  some  of  these  assumptions  is  to  use  a 
different  method  of  solving  the  acoustic  wave  equation. 

A  careful  analysis  of  the  errors  inherent  in  the  standard  PE 
method,  for  the  range- in depen dent  case,  is  carried  out  by  Fitzgerald 
C161.  The  importance  of  choosing  the  reference  sound  speed,  c^,  ,  to  be 
close  to  the  average  value  of  the  phase  velocities  c*,  for  excited 
modes,  is  quantified.  The  PE  errors  are  smaller  for  narrower  bands  in 
phase  velocities  or,  using  the  ray  mode  analogy,  the  narrower  the  cone 
of  propagating  rays.  Also,  it  is  shown  that  the  PE  method  maintains  a 
given  accuracy  to  longer  ranges  the  lower  the  frequency.  Since 
computation  time  increases  dramatical ly  with  increasing  frequency  in 
the  PE  method,  this  indicates  that  as  low  a  frequency  as  practical 
should  be  employed.  Also,  Fitzgerald  derives  a  compact  expression  for 
the  maximum  range  of  effectiveness  of  the  PE  method.  When  this  is 
applied  to  a  frequency  of  100  Hz,  in  a  range- independent  environment, 
with  a  simple  surface  and  bottom,  he  estimates  a  maximum  effective 
range  of  11  km.  For  a  10  Hz  source  located  near  the  surface,  so  that 
the  RSR  modes  are  the  only  excited  propagating  modes,  this  increases 
dramatically  to  15000  km.  Overall,  Fitzgerald  shows  that  the  PE 
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method  properly  accounts  for  mode  coupling  when  backscatter ing  is 
negligible,  and  properly  accounts  for  mode  propagation  when  there  is  a 
small  spread  in  mode  wavenumbers-  He  suggests  enhanced  PE  methods 
which  might  work  when  these  conditions  fail,  but  does  not  discuss  the 
greatly  increased  computation  time  that  would  likely  be  required  to 
carry  out  these  interesting  suggestions. 

A  different  approach  is  considered  by  DeSanto  C133,  who 
relates  the  solutions  of  the  general  Helmholtz  equation  and  the  PE  via 
an  integral  transform,  which  is  exact  for  range-independent  sound 
speeds.  However,  there  is  no  estimate  given  for  the  error  introduced 
when  the  Helmholtz  equation  is  replaced  by  the  PE.  It  appears  that 
DeSanto’s  interesting  ideas  could  lead  to  such  estimates,  which  would 
be  important  if  they  could  be  obtained. 

An  alternative  approach  to  the  numerical  computation  of  the 
wide-angle  acoustic  field  in  th©  range-dependent  case  is  given  in  C 13 1 
by  Estes  and  Fain.  They  also  proceed  stepwise,  by  employing  a  fully 
defined  field  in  one  vertical  plane  to  develop  the  field  at  a  plane  a 
fixed  distance  down  range.  The  propagation  is  composed  of  two  parts, 
somewhat  analagoue  to  the  predictor/cGrrector  methods  in  differential 
equations:  propagation  through  a  homogeneous  space,  and  a  correction 
due  to  the  fact  that  the  environment  is  not  homogeneous.  They 
approximate  the  square  root  operator  by  the  terms  of  its  Taylor  series 
up  to  order  4,  but  the  results  have  not  been  coded,  so  numerical 
evidence  is  unavailable.  Although  their  method  seems  to  avoid  the 
phase  errors  that  are  an  inherent  part  of  the  standard  (narrow-angle) 
PE  approach,  and  thus  offers  an  improvement  over  that  method. 
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appeals  that  the-  recent  development  of  the  wide-angle  PE  supercedes 
the  results  of  Estes  and  Fain.  Another  attempt,  the  first  order 
continued  fraction  approximation  of  Berkhout  C43,  is  equivalent  to  the 
standard  narrow  angle  approximation,  anc  his  second  order  continued 
fraction  result  is  just  the  Claerbout  (u/ide-angle)  approximation. 

Fitzgerald  C161  discusses  the  possibility  of  employing  normal 
mode  theory  to  handle  the  range-dependent  case,  by  dividing  the  range 
into  intervals  within  each  of  which  the  sound  speed  is  assumed  to  be 
range-independent.  Note  that,  in  view  of  the  range- in dependence  of 
the  PE  method  described  above,  this  technique  is  also  required  (and  is 
used  extensively)  when  applying  the  P£  method  to  the  range-dependent 
case.  This  piecewise  constant  normal  mode  approach  is  not  analyzed 
by  Fitzgerald.  Other  passible  approaches  to  solving  the  acoustic  wave 
equation  will  be  discussed  below, 

10.  Solution  Methods  for  the  Parabolic  Equation 

The  original  numerical  method  of  salving  the  PE  is  known  as 
the  split-step  Fourier  method.  The  basic  equation  of  the  split-step 
Fourier  method  is  an  exact  solution  to  the  PE  only  for  constant  index 
of  refraction  n.  Since  vertical  sound  speed  changes  are  usually  less 
than  5%,  and  strong  near-surface  horizontal  sound-speed  gradients  are 
roughly  0.1  o/s/km,  the  split-step  method  will  often  give  a  valid 
solution  to  the  wave  equation  for  a  real  environment. 

The  split-step  method  of  solving  the  PE  is  fast  and  accurate 
when  bottom  effects  are  weak,  but  less  efficient  when  bottom 
interaction  is  strong.  A  viable  alternative  is  the  implicit  finite 
difference  (IFD>  method,  which  is  usually  comparat i vely  fast.  and 
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unconditionally  stable.  IFD  is  applicable  when,  as  the  PE  solution  is 
marched  out  in  range,  the  boundary  information  at  the  advanced  range 
step  cart  be  expressed  in  terms  of  the  known  values  at  the  present 
range  step.  The  Lee-Botse<*s  IFD  model  C 25 3  is  applicable  to  both 
range-dependent  and  range-independent  environments,  and  can  handle 
arbitrary  surface  and  bottom  boundary  conditions,  horizontal 
interfaces  of  layered  media  (with  densities  constant  on  each  side  of 
the  interface;  see  C29  3  and  [313  for  details  of  the  treatment  of 
horizontal  interfaces),  attenuation,  shallow  or  deep  water,  and 
shallow  to  deep  or  deep  to  shallow  propagation.  As  a  matter  of 
necessity  (as  explained  above),  in  a  range-dependent  environment  the 
SVP  can  be  updated  only  abruptly  at  the  range  steps,  and  cannot  change 
within  a  stepr.  Advantages  of  finite-difference  methods  over  split- 
step  methods,  in  some  idealized  range- in depen dent  cases,  are  discussed 
by  McDaniel  C30  3.  A  complete  analysis  of  the  wel 1-posedness  of  the 
IFD  scheme  for  the  solution  of  the  wide-angle  PE  is  given  by  St.  Mary 
and  Lee  [403. 

The  split-step  algorithm  solves  a  pure  ini t ial -value  problem, 
by  imposing  an  artificial  zero  bottom  boundary  condition  and  pressure- 
release  surface  condition.  The  implicit  finite  difference  (IFD) 
scheme,  on  the  other  hand,  is  designed  to  treat  the  bottom-boundary 
condition  exactly.  Furthermore,  since  strong  bottom  interaction  is 
often  associated  with  a  wide  propagation  angle,  the  wide-angle  IFD 
scheme  C9;  could  prove  to  be  especially  useful  in  handling  cases  where 
both  of  these  phenomena  occur. 

The  method  of  lines  has  also  been  employed  to  solve  the  PE. 
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This  technique  is  described  by  Lee  and  Papadakis,  who  call  it  the  ODE 
method.  One  advantage  of  this  approach  is  that,  if  the  bottom 
boundary  conditions  are  known,  no  subbottom  knowledge  is  required. 
The  ODE  technique  needs  considerable  computer  memory  capacity,  but  if 
this  is  available  the  method  is  efficient.  The  computer  memory 
requirement,  in  fact,  can  be  reduced  by  applying  the  Generalized  Adams 
Method  (GAM)  described  by  Lee,  Jackson  and  Preiser  C271.  They  claim 
that  GAM  reduces  computation  speed  tremendously  as  well.  McDaniel, 
Saad  and  Lee  C32  3  also  discuss  the  GAM. 

V.  Bottom  Conditions 

The  broad  effect  of  bottom  loss  on  ocean  sound  propagation  is 
an  increasing  loss  with  decreasing  frequency.  In  actuality,  the 
bottom  is  constituted  of  many  sedimentary  layers,  each  acting 
acoustically  like  a  fluid.  The  sound  speed  in  the  bottom  is  usually 
close  to  that  of  the  water,  but  the  attenuation  is  much  greater  in  the 
bottom.  We  describe  some  of  the  ways  which  have  been  devised  to  treat 
the  various  types  of  bottom  conditions  which  can  occur. 

One  method  of  representing  the  bottom,  chosen  by  Guthrie  and 
Gordon  1183,  is  as  a  series  of  fluid  layers,  where  in  the  ith  layer 
the  sound  speed  varies  as 

c( z )-cj C l-2gj < z-z j ) /c; 3  , 

where  c*  and  g$  are  the  sound  speed  and  sound  speed  gradient  at  the 
top  of  the  ith  layer,  and  Zj  is  the  layer  depth.  Their  specific 
example  has  5  layers,  with  the  fifth  one  having  a  small  negative 
gradient.  This  is  necessary  to  satisfy  the  mechanics  of  the  normal 
mode  starter  program  of  the  PE  method,  and  has  no  effect  on  the 
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solution,  as  its  depth  is  below  that  to  which  significant  acoustic 
energy  penetrates.  A  variable  jump  discontinuity  in  the  sound  speed, 
of  11%  on  the  slope  and  tapered  to  zero  on  reaching  the  ocean  basins, 
is  applied  at  the  water-sediment  interface.  This  allows  for  greater 
bottom  reflectivity  an  the  continental  slope  as  compared  with  the 
ocean  basins,  and  compensates  far  the  inability  to  include  density  in 
the  model.  Finally,  they  choose  a  constant  attenuation  coefficient  in 
the  bottom  of  0. If  dB/km  <f=frequency  in  Hz),  corresponding  to  a 
porosity  of  at  least  65%. 

It  is  obvious  from  this  description  that  Guthrie  and  Gordon 
carried  out  a  very  detailed  analysis  of  the  bottom  structure  and 
properties  before  applying  the  PE  method.  As  the  PE  method  is  quite 
sensitive  to  errors  in  the  knowledge  of  the  bottom  environment,  a 
similar  bottom  analysis  will  be  required  for  any  application  of  the  PE 
method  to  propagation  loss  calculations  to  succeed. 

When  the  split-step  method  is  used  to  solve  the  PE,  the  bottom 
conditions  must  be  modeled  in  such  a  way  that  the  potential  function 
goes  smoothly  to  zero  at  maximum  depth.  One  way  to  accomplish  this  is 
to  extend  the  bottom  so  deep  that  the  sound  propagating  downwards  does 
not  reach  the  lower  boundary  within  the  horizontal  range  being 
considered.  This  is  not  practical,  however,  because  of  the  large 
increase  in  the  number  of  necessary  FFT  points  that  results.  Another 
possibility  is  to  introduce  a  sufficiently  high  bottom  absorption,  but 
it  can’t  be  arbitrarily  large,  since  it  must  represent  the  actual 
attenuation.  Mathematically  this  is  accomplished  by  incorporating  a 
volume  absorption  coefficient,  which  increases  rapidly  with  depth 


until  the  acoustic  field  becomes  negligible,  into  the  PE,  In  practice 
a  combination  of  these  two  methods  can  be  employed,  with  the  depth 
extension  and  the  absorption  chosen  to  fit  the  particular  environment 
being  considered.  Another  method  of  accounting  for  bottom  absorption 
is  to  add  a  small  imaginary  part  to  the  wave  number  in  the  bottom. 

In  many  ocean  environments  differences  in  temperature, 
pressure  and  salinity,  as  well  as  periodic  deposits  of  sedimentary 
material,  create  a  layered  medium  in  both  the  water  and  the  bottom. 
Various  attempts  have  been  made  to  adapt  the  PE  method  to  treat  the 
interfaces  between  such  layers.  In  several  papers  128,  29,  313, 
McDaniel  and  Lee  use  finite  difference  techniques  to  handle  such 
boundary  conditions.  One  drawback  of  this  work  is  that  the  mesh 
spacing  is  uniform,  introducing  computational  difficulty  when  there 
are  mesh  points  which  do  not  lie  exactly  on  the  interface  boundary. 
This  problem  was  addressed  by  deG  Gribble  C123,  who  employed  a 
variable  mesh  spacing.  An  implicit  finite  difference  (IFD) 
computational  technique,  incorporating  the  horizontal  interface  ideas 
of  McDaniel  and  Lee  and  the  wide  angle  equation  of  Claerbout  Ill  3,  is 
currently  undergoing  extensive  testing.  It  remains  unclear  how 
satisfactory  these  methods  are,  since  substantial  experimental 
evidence  is  lacking. 

Although  the  PE  method  does  not  take  the  density  change  at  the 
water-sediment  interface  into  consideration,  so  that  a  somewhat 
greater  bottom  reflection  loss  than  normal  is  to  be  expected,  this  is 
unlikely  to  be  a  source  of  serious  error.  For  example,  a  density 
contrast  of  1.25:1  causes  a  decrease  in  reflection  loss  at  low  grazing 
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angles  of  less  than  0.5  dB/bounce,  while  a  57.  jump  in  sound  speed  at 
the  sediment  interface  causes  a  more  meaningful  A  to  5  dE  decrease 
C  183. 


VI.  Ice  and  Other  Environmental  Concerns 

The  importance  of  accurate  knowledge  of  environmental 
conditions,  when  doing  propagation  loss  calculations,  is  abundantly 
clear  in  the  results  described  by  Ellis  and  Chapman  C141  for  low- 
frequency  shallow  water  propagation.  In  one  experiment,  curves  of 
propagation  loss  versus  frequency  at  a  range  of  55  km  are  given  for 
two  sites;  one  with  a  chalk  bottom  and  the  other  a  nearby  site  with  a 
sand  bottom.  For  frequencies  less  than  100  Hz,  the  loss  observed  over 
the  chalk  bottom  was  more  than  *5  dB  greater  than  that  over  the  sand 
bottom!  Thus,  certainly  for  shallow  water  low-frequency  propagation, 
knowledge  of  the  bottom  conditions  is  absolutely  essential  if  accurate 
calculations  are  to  be  made. 

A  triangular  ridge  model  has  been  shown  by  Greene  and  Stokes 
C173  to  give  better  agreement  with  backscatter  data  from  Arctic  sea 
ice  than  a  standard  Gaussian  model.  Their  model  and  data  confirm  the 
observations  of  Mellen  and  Marsh  C331  and  Brown  and  Milne  C7,  8,  341, 
that  such  backscatter  levels  are  much  stronger  than  those  in  the  open 
ocean.  This  is  caused  by  the  large  linear  ridge-keel  structures  which 
interrupt  the  otherwise  smooth  stretches  of  ice  on  the  underside 
surface  of  Arctic  sea  ice.  Furthermore,  small-scale  roughness,  due  to 
the  rubble  which  makes  up  the  keels,  is  superimposed  on  the  linear 
structure. 
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Since  a  fundamental  premise  in  the  PE  method  is  that  there  is 
no  backscattering,  it  is  clear  from  the  above  that  propagation  loss 
calculations  employing  the  PE  method  will  likely  contain  very  large 
errors  in  an  Arctic  underice  environment.  One  possible  mitigating 
factor  could  be  the  use  of  low  frequencies,  as  the  local 
backscatter ing  strength  of  a  patch  on  a  ridge  is,  according  to  Bass 
and  Fuks  C2,  p.  1081,  proportional  to  the  fourth  power  of  the  acoustic 
frequency.  Since  mc^t  of  the  published  data  involves  frequencies 
greater  than  1  kHz  (Greene  and  Stokes  C17  3,  for  example,  exhibit  plots 
of  backscattering  strength  at  1.81  kHz),  the  possibility  remains  that 
the  PE  method  will  give  reasonably  accurate  results  in  the  10B  Hz 
range.  Another  factor  to  consider  in  this  regard  is  that  PE 
approximations  are  better  at  low  grazing  angles,  which  is  precisely 
the  case  for  which  the  backscatter ing  probability  is  greatest.  To 
more  fully  answer  these  questions,  it  will  be  necessary  to  carry  out 
detailed  and  precise  Arctic  experiments,  similar  to  those  undertaken 
in  the  South  Tasman  Sea  in  1975  and  described  by  Guthrie  and  Gordon 
C  183. 

The  knowledge  of  bottom  properties  and  gross  sound-speed 
character istics  alone  will  probably  not  be  sufficient  to  obtain  an 
accurate  model  of  acoustic  propagation  under  ice.  For  example,  there 
may  be  energy  loss  as  a  result  of  mode  coupling,  either  from  ice- 
bottom  roughness  where  there  is  no  deep  sound  channel,  or  from 
inhomogeneities  within  the  water  column  in  the  boundary  region  between 
different  water  mass  types.  The  latter  could  be  caused,  for  example, 
by  marked  changes  in  salinity  as  the  sound  propagates  past  partially 


melted  ice  flows. 


Since  rays  which  propagate  at  steep  angles  are  greatly 
attenuated  by  radiation  and  bottom  absorption,  low  frequency  long 
range  propagation  is  dominated  by  rays  having  small  grazing  angles. 
In  this  environment,  therefore,  the  standard  PE  is  usually  acceptable. 
However,  when  the  reference  wave  number  differs  appreciably  from  the 
wave  numbers  of  the  propagating  modes,  phase  errors  are  introduced  in 
the  standard  PE.  These  errors  are  not  reduced  by  reducing  the  range 
step  or  increasing  the  number  of  FFT  points  in  the  split-step  solution 
of  the  PE.  A  major  advantage  of  the  wide  angle  PE  is  that  these  phase 
errors  are  essentially  eliminated. 

The  PE  method  has  been  applied  to  wide  angle  propagation  by 
employing  a  more  accurate,  rational  approximation  to  the  square  root 
operator  Q  defined  above.  Claerbout  Ill]  introduced  this  means  of 
attack,  by  using  ( l+3x/4)/(  l+x/4)  to  approximate  (1+xA,  in  place  of 
the  estimate  l+x/2  which  results  in  the  standard  narrow  angle  PE.  In 
this  context,  "narrow  angle"  means  less  than  about  15®,  while  "wide 
angle*  is  less  than  about  40®.  A  "very  wide  angle"  (less  than  about 
70*)  PE,  which  utilizes  a  second  order  rational  approximation  for  the 
square  root  operator,  also  exists,  but  its  utility  has  not  been 
demonstrated. 

Boundary  fluctuations,  source  directional ity  and  water  column 
variations  are  three  mechanisms  which  can  introduce  azimuthal 
variation  in  the  sound  pressure  field.  The  development  of  three- 
dimensional  (3D)  capability  in  the  various  PE  models  has  only  recently 
been  begun,  and  the  results  of  realistic  test  examples  are  not 


available  yet.  Siegmann  and  Lee  C  39  1  offer  a  more  detailed 
explanation  of  why  such  techniques  are  important.  They  also  develop 
estimates  of  ranges  for  which  a  3D  PE  would  be  valid.  For  example, 
their  estimates  mean  that  at  150  Hz  the  range  must  be  at  least  40  m, 
while  it  must  be  at  least  400  m  when  the  frequency  is  15  Hz. 

VII.  Choice  of  Parameters 

Uhen  studying  ducted  propagation  in  the  ocean,  there  is 
generally  an  optimum  frequency  to  consider.  Such  a  phenomenon  occurs 
because  of  competing  attenuation  and  propagation  mechanisms  at  various 
frequencies.  At  high  frequencies,  scattering  and  volume  loss  simply 
increase  with  frequency,  but  at  low  frequencies  the  situation  is  more 
complex.  Here  tha  intensity  of  sound  within  the  duct  is  affected  by 
propagation  and  attenuation  mechanisms  in  the  bottom  and  otherwise 
outside  the  duct.  As  sound  penetration  into  the  bottom  increases  with 
decreasing  frequency,  the  overall  attenuation  of  sound  throughout  the 
water  column  will  increase  as  the  frequency  decreases.  Therefore  both 
high  and  low  frequencies  can  result  in  high  attenuation,  with 
intermediate  frequencies  having  the  lowest  attenuation.  For  example, 
an  optimum  frequency  on  the  order  of  hundreds  of  hertz  occurs  in  a 
typical  shallow  water  environment. 

The  importance  of  the  choice  of  receiver  depth  is  made  clear 
in  the  results  of  Jensen  C201.  He  discusses  a  simple  example  of 
propagation  in  shallow  (110  m)  isothermal  water,  modeled  by  normal 
mode  theory,  at  a  frequency  of  100  Hz.  In  this  case  tha  theory  and 
experimental  data  agree  extremely  well.  Both  show  a  propagation  loss 
level  that  is  20-25  dB  greater  for  a  receiver  depth  of  1  m  than  that 
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for  a  receiver  depth  of  50  m,  over  the  entire  range  of  30  km. 

The  ideal  vertical  step  size  in  applying  the  split-step 
Fourier  algorithm  to  numerically  solve  the  PE  is  0.2^sint>-,  inhere  J*-  is 
the  wavelength  and  &  is  the  grazing  angle.  This  places  an  upper  limit 
on  the  number  of  modes,  and  hence  the  frequency,  that  can  be  handled, 
since  computer  run-time  is  proportional  to  n  log(n),  where  n  is  the 
number  of  FFT  points.  A  deterministic  criteria  for  the  ideal 
horizontal  step  size,  ir,  does  not  appear  to  be  known,  although  it  is 
certainly  limited  by  the  frequency,  and  the  magnitude  of  horizontal 
gradients.  In  practice  a  trial  and  error  approach  is  usually  taken, 
where  6r  is  decreased  until  the  solution  converges. 

It  is  essential  to  the  PE  method  that  the  initial  pressure 
distribution  over  depth  includes  full  phase  information  (coherent 
addition  of  modes).  Three  methods  of  starting  the  PE  model  are  normal 
mode  theory,  ray  theory,  and  assuming  a  Gauss ian-shaped  initial 
pressure  distribution  near  the  source.  The  latter  is  appropriate  for 
deep  water,  because  it  filters  out  high-angle  energy  (which  would  be 
lost  anyway)  and  prevents  aliasing  of  this  energy  into  low  angles 
during  the  FFT  sampling  process.  Initializing  the  PE  model  by  normal 
mode  theory  requires  a  constant  water  depth  for  the  interval  in  which 
the  mode  theory  is  being  carried  out,  and  imposes  an  upper  limit  of 
about  60  degrees  on  the  equivalent  ray  angle.  The  ray  starter  seems 
best  for  propagation  down  the  continental  slope  and  into  deep  water. 
In  this  case  very  high  angle  energy  close  to  the  source  must  be 
considered,  because  of  the  considerable  reduction  in  the  angle  during 
propagation  down  the  slope.  On  the  basis  of  the  requirement  that 
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energy  which  insanities  convergence  zones  when  it  reaches  deep  water 
must  be  treated  properly  at  the  top  of  the  slope,  Guthrie  and  Gordon 
C 18 1  choose  a  vertical  sample  interval  of  Jz=2H/5m,  where  H  is  the 
water  depth  and  m  is  the  highest  order  mode  required  in  deep  water. 
Since  at  100  Hz  it  is  necessary  to  consider  about  150  modes  in  deep 
water,  a  step  size  of  12  m  is  called  for  in  4500  m  water.  At  the  top 
of  the  continental  slope,  where  H  is  about  350  m,  a  step  size  of  about 
1  m  is  adequate. 

In  the  extensive  comparisons  of  experimental  data  with  PE 
model  predictions  described  by  Guthrie  and  Gordon  C183,  at  frequencies 
of  63  and  125  Hz,  the  range  step  size  Jr  varied  from  10  m  on  the 
continental  shelf  to  200  m  in  deep  water,  while  the  depth  step  size  Jz 
varied  from  6  m  to  12  m,  with  the  greater  Az  also  taken  in  deep  water. 
However,  they  speculated  that  under  certain  circumstances  a  Az  of  3  m 
or  less  would  have  improved  the  results,  but  such  an  option  was  not 
available  to  them  in  their  model.  This  indicates  that  the  choice  of 
vertical  step  size  must  be  made  with  considerable  care.  Az  must  be 
small  enough  to  ensure  that  the  highest  order  modes  needed  are  sampled 
properly,  but  not  so  small  that  computer  run-time  becomes  excessive. 
In  the  interest  of  efficiency,  Jz  should  be  chosen  automatical ly  by 
the  model  to  suit  local  conditions. 

Jensen  and  Krol  C19J  carry  out  a  detailed  analysis  of  the 
allowable  range  step,  Jr,  when  the  split-step  Fourier  method  is  used 
to  solve  the  PE.  They  conclude  that  for  deep  water  Jr  can  be  taken  as 
large  as  100  or  even  1000  m.  In  general,  it  is  required  to  take  Jr 
inversely  proportional  to  both  the  maximum 
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gradient  and  to  the  frequency*  Furthermore,  is  slightly  dependent 
upon  bottom  attenuation  and  source/receiver  depths.  They  recommend, 
as  did  Guthrie  and  Gordon,  that  4r  for  a  given  propagation  problem  be 
determined  empirically,  by  reducing  it  until  a  stable  solution  has 
been  obtained.  Finally,  Jensen  and  Krol  observe  that  in  a  shallow 
water  environment  it  appears  that  far  must  be  extremely  small. 

VIII.  Comparisons  With  Other  Models  and  Experimental  Data 

In  order  to  determine  the  accuracy  and  applicability  of  the  PE 
method,  it  must  be  compared  with  both  experimental  data  and  with  the 
predictions  resulting  from  other  models.  While  agreement  with  other 
models  instills  confidence  in  any  particular  method,  the  final  check 
on  an  acoustic  model  is  a  comparison  with  the  results  of  realistic 
experiments.  Only  then  should  we  be  convinced  that  the  model  being 
evaluated  includes  all  of  the  necessary  and  correct  physics  and 
mathematics  for  understanding  and  explaining  sound  propagation  in  a 
real  ocean  environment. 

An  example  of  a  complicated  North  Atlantic  environment  modeled 
by  adiabatic  normal  mode  theory  is  presented  by  Jensen  C203.  The 
water  depth  varies  between  113  and  305  m,  with  a  maximum  bottom  slope 
of  about  1*.  The  65  km  range  is  divided  into  two  distinctly  different 
water  masses,  with  two  different  SVP’s.  The  bottom  is  also  of  two 
types,  sand  and  silt,  with  the  changeover  occurring  at  a  range 
different  from  that  where  the  SVP  changes.  In  spite  of  the  difficult 
environment,  Jensen  found  excellent  agreement  between  theory  and 
experiment.  Over  almost  6  octaves  <50-2540  Hz)  of  frequency,  and  over 
the  entire  range  of  65  km,  the  maximum  deviation  was  only  a  few 
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If  the  calculation  times  (which  are  not  mentioned)  required  to  obtain 
these  results  are  at  all  reasonable,  this  strongly  suggests  that 
adiabatic  normal  mode  theory  is  a  very  viable  alternative  to  the  PE 
method  for  a  range— dependent  environment.  It  should  be  emphasised, 
however,  that  in  this  example  an  extremely  detailed  environmental 
description  played  a  crucial  role  in  the  achieved  accuracy  of  the 
modelling  results. 

A  detailed  experiment  of  propagation  over  a  seamount  is 
described  by  Jensen  (C  20  3,1211),  with  the  PE  method  used  as  a  model. 
Although  he  calls  this  a  'strongly  range-aependent  environment  <14® 
mean  slope  across  the  mount),"  the  SVP  remains  constant  over  the 
entire  140  km  range.  Thus  the  common  terminology  would  refer  to  this 
as  a  range— independent  case.  At  any  rate,  the  agreement  between 
theory  and  experiment  is  uniformly  good  over  the  full  14©  km  range, 
for  frequencies  ranging  from  12.5  to  408  H2.  Moreover,  in  this  case 
the  agreement  is  achieved  without  very  detailed  environmental 
information.  This  indicates  that  the  real  SVP  is  well  approximated 
by  the  constant  one  that  is  assumed,  and  that  the  important  physics  is 
associated  with  reflection  off  the  seamount  and  with  propagation  in 
the  water  column.  Thus  these  results  can  be  somewhat  encouraging  for 
the  study  of  under-ice  propagation,  although  backscatter i ng  off  ice 
keels  and  SVP  changes  due  to  salinity  differences  are  still  likely  to 
present  significant  difficulties,  as  described  earlier. 

Comparisons  of  normal  mode,  fast  field  program  (FFP),  and  PE 
wave  theory  models  for  test  cases  involving  four  different  underwater 
acoustic  environments  are  described  in  1223  by  Jensen  and  Kuperraan.  A 
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split-step  algorithm  (PAREQ)  is  used  to  solve  the  PE-  The  first  case 
involves  a  range-dependent  environment,  where  the  SVP  remains  constant 
for  the  first  20  kro,  transitions  for  10  km,  and  remains  constant  again 
for  the  final  20  km.  The  frequency  is  25  Hz.  For  receiver  depths  of 
25,  250,  and  400  m  the  three  models  agree  quite  well  over  the  first  20 
km.  The  normal  mode  solution  is  stopped  at  this  point,  because  of  the 
range-dependent  SVP.  However,  when  the  FFP  and  PE  solutions  are 
continued  out  to  the  maximum  range  of  50  km,  they  diverge 
substantially  from  each  other.  Jensen  and  Kuperman  argue  for 
accepting  the  PE  solution,  stating  that  "the  PE  gives  the  correct 
result  in  the  range-independent  part  of  the  environment  and  because  of 
the  mathematical  nature  of  the  PE  solution  there  is  no  reason  to 
believe  that  the  accuracy  of  this  solution  diminishes  in  the  mildly 
range— dependent  region.*  We  question  this  argument,  because  of  the 
proof  given  earlier  that  the  derivation  of  the  PE  implies  a  range- 
independent  SVP. 

The  second  case  in  C  22  3  is  a  range- independent  environment 
with  different  bottom  speeds.  The  frequency  is  again  25  H2.  As*  the 
maximum  propagation  angle  increases  with  increasing  bottom  speed,  and 
as  the  standard  P€  is  valid  only  for  narrow  angles  of  propagation,  the 
agreement  between  the  PE  solution  on  the  one  hand,  and  the  apparently 
accurate  FFP  and  normal  mode  solutions  on  the  other  hand,  deteriorates 
with  increasing  bottom  speed.  This  offers  a  good  illustration  of  the 
need  for  the  wide  angle  approximation  to  the  PE,  discussed  elsewhere 
in  this  paper. 

Case  3  in  C  22  3  is  a  range- independent  shallow  water  (100  m) 
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environment,  at  a  frequency  of  250  Hz.  With  the  source  and  receiver 
depths  both  50  m,  the  agreement  between  the  three  models  is  excellent 
over  the  full  range  of  10  km,  except  for  3  unexplained  isolated  spikes 
in  the  PE  curve,  which  overpredict  the  propagation  loss  by  as  much  as 
35  dB.  When  the  source  and  receiver  are  located  just  off  the  bottom, 
so  that  higher  order  modes  are  excited,  the  PE  solution  no  longer 
agrees  as  well  with  the  correct  normal  mode  and  FFP  solutions.  The 
reason  for  this  is  the  increased  inaccuracy  in  the  reference  sound 
speed  when  more  snodss  are  excited,  as  we  discuss  elsewhere. 

Ths*  fourth  case  in  Z  22 1  deals  with  a  fairly  complicated  range- 
independent  environiaent,  where  the  bottom  is  sloping  and  has  a 
sediment  layer.  The  frequency  is  25  Hz.  Agreement  between  the  three 
models  is  quite  good,  even  out  to  a  range  of  150  km. 

In  addition  to  the  results  of  the  four  test  cases  Just 
described,  Jensen  and  Kuperman  1223  list  the  basic  numerical 
parameters  and  associated  computer  times  for  the  normal  mode  <NH)  and 
PE  solutions.  The  depth  step  varies  from  0.5  m  (case  3)  to  0.5  m 
(case  4 >  for  NM,  and  from  0.26  m  (case  3)  to  11.2  >  (case  2)  for  PE. 
The  range  step  4r  varies  from  10  m  (case  3)  to  500  m  (case  4)  for  NM, 
and  from  2  m  (case  3)  to  100  m  (cases  1  and  4)  for  PE.  The  number  of 
FFT  points  for  the  PE  is  256  for  case  2,  and  1024  for  cases  1,  3  and 
4.  The  calculation  time  on  a  UNIVAC  1100/60  varies  from  ’ess  i han  1 
min  (cases  2  and  3)  to  2  min  (case  4)  for  NM,  and  from  1  min  (case  2) 
to  30  min  (case  3)  for  PE.  The  PE  method  is  particularly  slow  when 
there  is  heavy  bottom  interaction  (30  min  in  case  3  and  5  min  in  case 
4) . 
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Examples  of  the  application  of  the  IFD  technique  to  several 
problems  are  described  by  Lee  and  Botseas  C25  3,  and  compared  to 
solutions  obtained  by  normal  mode  theory,  and  by  Jensen  and  Kuperman 
C223  using  the  split-step  method  of  solving  the  PE.  In  general  the 
solutions  agree  quite  well.  The  reader  is  referred  to  C253  for  the 
details. 

Guthrie  and  Gordon  t  18  3  describe  extensive  propagation  loss 
calculations  and  experiments  carried  out  in  the  south  Tasman  Sea  in 
1975.  This  appears  to  be  the  most  definitive  3tudy  to  date  on  the 
appropriateness  of  the  PE  method  in  a  realistic  ocean  environment. 
The  computer  runs  involve  accurate  bathymetric  and  sound-speed  data 
for  ranges  up  to  2503  km  at  frequencies  of  63  and  125  Hz.  Among  the 
interesting  and  important  results  obtained  were: 

J.  Initialization  of  the  PE  model  in  a  shallow,  coastal 
environment  was  accomplished  by  a  ray  trace  and  a  normal  mode 
prediction,  and  in  both  cases  the  models  predicted  absolute  levels  in 
deep  water  that  were  within  5  dS  of  those  obtained  experimental ly. 

2.  The  optimum  vertical  sample  size  was  difficult  to 
determine  when  running  a  P£  prediction  from  shallow  to  deep  water,  or 
vice  versa. 

3.  The  assumed  properties  of  the  sediment  cover  had  a  marked 
effect  upon  the  predicted  acoustic  shadow  of  substantial  underwater 
features. 

9.  The  model  failed  to  predict  an  observed  increase  in  shadow 
depth  with  frequency. 

5.  The  clarity  of  the  convergence  zone  structure  varied 
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strongly  with  range,  in  certain  range-dependent  environments. 

6.  The  predicted  acoustics  were  very  strongly  dependent  upon 
the  SVP,  with  even  subtle  changes  in  sound-speed  profiles  having  a 
marked  effect. 

7.  At  very  low  frequencies  <16  and  31.5  Hz),  the  relatively 
weak  source  level  and  high  ambient  noise  conditions  resulted  in  poor 
quality  data. 

Several  other  detailed  comparisons  of  experimental  data  and  PC 
predictions  are  given  by  Guthrie  and  Gordon  [18  3.  The  reader  is 
referred  to  the  original  paper  (which  should  be  required  reading  for 
anyone  in  the  field)  for  the  precise  results.  W©  only  observe  here 
that,  although  agreements  between  experiment  and  prediction  are 
usually  reasonably  good,  discrepancies  on  the  order  of  2.0  dB  or  more 
are  quite  common.  Moreover,  on  occasion  the  errors  were  so  large  as 
to  render  the  particular  PE  model  being  employed  useless,  and 
requiring  ad  hoc  adjustments  to  the  model  (such  as  changes  in  the 
bottom  assumptions)  to  bring  the  predictions  buck  into  line  with  the 
data.  This  clearly  indicates  the  importance  of  achieving  accurate 
knowledge  of  the  propagation  character ist ics  of  the  region  in  which 
one  wants  to  employ  the  PC  method,  before  the  actual  work  is  carried 
out. 

Jensen  and  Krol  C193  give  propagation  loss  predictions  for 
several  cases,  employing  the  PE  split-step  model,  a  normal  mode  model, 
and  a  ray-tracing  model.  Por  a  deep  water  (3000  m)  range- in depen dent 
case,  with  a  frequency  of  100  Hz  and  a  range  of  65  km,  the  curves 
generated  by  the  three  models  agree  very  well.  The  range  step  dr  for 
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the  PE  program  is  15300  m.  On  a  UNI VAC  1106,  the  computation  times 
were  6  min  for  the  ray  program,  12  min  for  the  mode  program,  and  24 
min  for  the  PE  program.  For  a  deep  water  range-dependent  case,  with  a 
frequency  of  250  Hz  and  a  range  of  32  km,  ray-tracing  and  PE 
predictions  are  compared  with  experimental  data,  and  in  both  cases 
agreement  is  quite  good.  4r  for  the  PE  program  is  100  m,  and 
calculation  times  are  4  hours  for  the  ray-tracing  program  and  28  min 
for  thr>  program.  For  a  shallow  water  range— independent  case,  the 
results  of  a  normal  mode  model  and  the  PE  model  agree  very  well. 
However,  a  4r  of  2  ®  is  required  to  obtain  a  stable  solution  for  the 
PE  model.  This  small  range  step  results  in  a  calculation  time  of  3.5 
hours  for  the  PE  program,  whereas  the  mode  program  takes  less  than  one 
min  to  execute. 

The  results  of  several  numerical  computations  are  shown  in 
C 1 3.  They  all  deal  with  the  case  where  the  SVP  is  range-independent, 
and  the  normal  mode  solution  is  used  as  a  reference.  For  a  narrow 
angle  (5.5*)  test  case,  agreement  between  the  normal  mode  solution  and 
the  PE  solution  employing  the  IFD  scheme  with  interface  is  quite  good. 
The  phase  error  in  the  standard  PE  is  evident  in  a  wider  angle  (19*) 
test  case,  but  the  wide  angle  PE  gives  excellent  agreement  with  the 
normal  mode  reference  solution.  For  a  very  wide  propagation  angle 
(70*),  no  available  PE  method  produces  a  satisfactory  result. 

Bates  C33  implemented  a  range-independent,  boundary-free 
example  with  an  SVP  chosen  so  that  the  wave  numbers  and  modal  depth 
functions  could  be  obtained  in  closed  form.  The  split-step  algorithm 
for  solving  the  standard  PE  was  then  tested,  by  comparison  with  the 
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closed  form  solution,  for  various  choices  of  frequencies,  number  of 
FFT  points,  and  sound  speed  gradients.  To  indicate  the  frequency 
dependence,  for  a  single  propagating  mode  at  50  Hz,  a  range  increment 
of  400  meters  and  64  FFT  points  mas  sufficient  to  generate  results 
uiithin  1  dB  of  the  closed  form  solution.  When  the  frequency  mas 
increased  to  500  Hz  (and  all  other  parameters  remained  the  same),  the 
error  increased  to  5  dB,  and  at  1  kHz  it  mas  11  dB.  When  the  number 
of  FFT  points  mas  varied,  with  the  frequency  held  constant  at  1  kHz 
and  1  propagating  mode,  &  doubling  in  the  number  of  FFT  points 
resulted  in  a  reduction  in  the  error  by  a  factor  of  4.  For  example, 
the  error  was  42  dB  with  32  FFT  paints,  11  dB  with  64  FFT  points,  2.75 
dB  with  120  FFT  points,  and  .67  dB  with  256  FFT  points.  As  expected, 
an  increase  in  the  sound  speed  gradient  resulted  in  an  increase  in  the 
error,  and  an  increase  in  the  number  of  propagating  modes  from  1  to  3 
caused  the  introduction  of  phase  errors.  It  is  interesting  to  note 
that  the  magnitude  of  these  phase  errors  was  independent  of  the  number 
of  FFT  points,  at  least  for  the  limited  number  of  examples  considered 
by  Bates.  Also  interesting  is  the  fact  that  all  of  the  errors 
discussed  here  were  essentially  independent  of  range,  for  those  ranges 
considered  (up  to  10  km). 

Results  of  a  test  of  the  wide-angle  1FD  model  are  reported  in 
C53.  The  test  case  is  a  range- in depen dent  environment,  with  an 
isovelocity  water  column  over  an  isovelocity  half-space  bottom.  The 
frequency  is  250  Hz,  the  water  depth  is  100  m,  the  maximum  range  is  IQ 
km,  and  the  source  and  receiver  are  just  off  the  bottom.  Because  of 
this  last  condition  the  higher  modes  are  more  strongly  excited,  anc 


all  11  modes  of  this  example  are  required  for  an  accurate  solution. 
The  reference  solution  is  one  gotten  from  the  fast-field  program 
(FFP)f  which  has  been  compared  with  normal  mode  results  and  is 
believed  to  be  correct.  In  this  test  the  wide-angle  IFD  results  are 
in  excellent  agreement  with  the  FFP  reference  solution.  When  the 
narrow-angle  (standard)  IFD  scheme  is  used,  the  expected  phase  errors 
appear . 

Comparisons  of  experimental  propagation  loss  data  with 
predictions  made  by  the  PE  method  and  by  ray  tracing  are  described  in 
161,  for  convergence  zone  propagation  over  60  km,  at  a  frequency  of 
3S0  Hz,  with  a  range-independent  SVP.  The  PE  solution  agreed  with  the 
measured  data  to  within  1  dB,  and  the  ray  trace  solution  agreed  within 
1.3  dB.  Note  that  in  this  case,  of  the  first  convergence  zone,  there 
is  very  little  bottom  interaction,  so  that  the  bottom  handling 
capabilities  of  the  PE  method  were  not  tested  in  this  experiment. 
Furthermore,  when  multiple  SVP’s  measured  along  the  propagation  track 
were  input  to  the  PE,  and  the  resulting  predicted  propagation  loss  was 
compared  with  experimental  data,  the  comparison  was  not  significantly 
better  than  that  obtained  previously  by  employing  a  single  SVP  over 
the  entire  track.  Thus,  incorporating  SVP  range— dependence  did  not 
appreciably  improve  the  results. 

Chan  *t  al  £91  point  out  the  existence  of  an  unconditionally 
stable  IFD  scheme,  which  discretizes  the  standard  3D  PE  by  means  of 
central  finite  differences  for  both  the  depth  and  azimuthal 
derivatives,  and  then  applies  the  Crank-Nicolson  method.  The  results 
appear  to  be  reasonably  accurate,  but  the  scheme  requires  excessive 
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computer  time.  In  C93  they  introduce  a  stability-generating 
dissipative  term  into  a  basic  finite  difference  scheme  for  the 
standard  Schrodinger  equation,  and  analyze  and  test  the  result.  In 
their  test  case  this  explicit  scheme,  using  the  same  range  step 
(0.081m)  as  in  the  Crank-Nicolson  implicit  scheme,  produces 
essentially  the  same  results,  but  3.5  times  faster  in  terms  of  CPU 
time.  The  actual  computation  times  are  not  mentioned.  The  accuracy 
of  the  computed  solutions  is  measured  by  comparing  them  to  what  is 
claimed  to  be  an  exact  solution,  u=sin(<<z)exp(  imO)exp< imV(2k*r ) ) . 
However,  it  is  straightforward  to  show  that  if  this  function  u 
satisfies  the  standard  3D  PE,  then v(n*-l ).  Since  <<  is  constant 
(TT/100  in  their  teat  case),  this  requires  that  n  is  also  constant, 
i.e.,  the  range-independent  case.  m  tact,  the  last  equation  gives  an 
explicit  relationship  between <*,  k* ,  and  n,  which  must  hold  for  this 
test  case  to  be  valid.  The  test  frequency  is  not  mentioned  in  C9  3,  so 
it  is  not  clear  if  this  consideration  was  taken  into  account.  The 
test  sound  speed,  however,  was  taken  to  be  constant,  by  necessity  (as 
explained  above),  not  for  Simplicity*  (as  claimed  in  C93). 

In  E103  Chan  et  al  describe,  in  somewhat  more  detail  than  in 
C9  3,  various  finite  difference  schemes  for  solving  the  P€.  They 
define  a  "practical  stability  condition,"  and  compare  thr*e  stable 
explicit  finite  difference  schemes  with  the  Crank-Nicholson  method. 
The  greatly  oversimplified  equation  ur=»(  i/2kj,  )u**  is  used  as  a  test 
case,  with  a  depth  of  160  m,  a  frequency  of  500  Hz,  a  constant  sound 
speed  of  1580  ta/s,  a  range  of  200  to,  a  depth  increment  of  1  m,  and  a 
range  step  of  0.001  m. 
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Apparently  this  small  range  step  was  required  in  order  for  the 
four  methods  to  agree  with  the  exact  solution  to  within  IX.  In  spite 
of  the  simple  equation  being  solved,  the  lack  of  variation  in  the 
environmental  parameters,  and  the  short  range,  the  CPU  times  required 
were  substantial.  For  the  three  explicit  schemes,  these  times  ranged 
from  more  than  48  minutes  to  almost  1  hour  and  22  minutes,  while  the 
Crank-Nicholson  method  took  over  2  hours  and  7  minutes  of  CPU  time. 

IX.  Possible  Advantages  of  the  Parabolic  Equation  Method 

Although  ray  theory,  particularly  at  higher  frequencies, 
generally  requires  much  less  calculation  time  than  wave  models  such  as 
the  PE  method,  ray  theory  is  incapable  of  handling  diffraction.  This 
causes  a  ray  theory  model  to  predict  zero  sound  in  shadow  regions 
which,  although  unrealistic,  might  be  satisfactory  for  the  many 
applications  where  the  primary  interest  is  in  sound  levels  in  the 
insonified  regions.  In  other  cases,  where  diffraction  is  important, 
the  PE  method  or  some  other  wave  model  must  be  used,  because  only  such 
models  accurately  predict  the  frequency  dependence  of  the  diffracted 
field  and  hence  give  correct  field  levels  in  the  shadow  zones. 

In  a  study  of  propagation  in  a  wedge-shaped  ocean  with  a 
penetrable  bottom  £231,  Jensen  and  Kuperman  observe  that  the  PE  method 
performs  a  reasonably  good  Job  of  modeling  the  modal  cutoff  that 
occurs  during  ups lope  propagation.  Adiabatic  mode  theory,  on  the 
other  hand«  fails  to  give  good  results  in  this  case,  because  of  the 
penetrability  of  the  bottom. 

In  the  PE  model  both  the  SVP  and  the  bottom  depth  and 
structure  are  allowed  to  vary  in  range.  Attenuation  in  the  bottom  is 


included,  by  introducing  a  complex  sound  speed. 

Although  the  effect  of  ocean  current  profiles  on  sound 
propagation  i®  at  least  an  order  of  magnitude  less  that  the  effect  of 
sound-speed  profiles,  the  cumulative  effect  of  current  profiles  may 
cause  significant  changes  in  PE  predictions  of  underwater  acoustic 
fields.  The  effect  of  current  on  reciprocity  is  more  pronounced.  It 
is  easy  to  numerically  include  current  profiles  in  the  PE  model,  with 
a  negligible  increase  in  program  complexity  or  calculation  time.  In 
cases  where  one  is  interested  in  absolute  levels  of  propagation  loss, 
the  ability  to  include  currents  is  certainly  an  important  feature. 
When  relative  las*  levels  are  of  concern  currents  would  probably  not 
matter  as  much,  since  they  would  have  the  same  effect  on  different 
sianals  of  the  same  freauency.  However,  current  considerations  might 
very  well  play  a  role  in  the  choice  of  search  depth  and  pitch  angle. 
Phu  and  Tappert  C333  give  a  detailed  discussion  of  why  and  how  current 
should  be  incorporated  in  PE  models. 

The  phase  errors  inherent  in  the  application  of  the  standard 
narrow-angle  PE  Vtave  been  essentially  eliminated  by  development  of  the 
wide-angle  PE,  The  key  to  the  derivation  of  the  wide-angle  PE  is  the 
approximation  of  the  square  root  operator,  which  naturally  arises  when 
the  PE  is  being  formulated,  by  a  nonlinear,  rational  function,  as 
opposed  to  the  linear  approximation  which  yields  the  standard  PE.  A 
detailed  derivation  of  the  wide-angle  PE  is  given  by  Siagraann, 
Kriegsmann  and  Lee  C  393.  They  include  an  error  analysis,  which  we 
believe  is  flawed,  as  described  in  this  paper.  A  split-step  algorithm 
for  the  wide-angle  PE  appears  in  the  paper  of  Thomson  and  Chapman 


[433.  They  point  out  that  existing  computer  codes  for  the  split-step 
solution  of  the  standard  PE  can  be  easily  modified  to  incorporate  the 
wide-angle  capability,  with  very  little  change  in  computation  time. 

The  PE  method  allows  the  acoustic  wave  equation  to  be  treated 
as  an  initial  value  problem,  which  is  amenable  to  a  marching  solution. 
In  contrast  the  elliptic  equation,  from  which  the  PE  is  derived, 
requires  boundary  conditions  which  must  be  prescribed  on  a  closed 
surface.  Thus  the  PE  method  is  computationally  much  more  efficient 
than  trying  to  numerically  solve  the  elliptic  acoustic  wave  equation. 
Also,  it  is  claimed  by  Tappert  and  Hardin  C423  that  the  split-step 
Fourier  method  of  solving  the  PE  is  unconditionally  stable,  even  when 
the  index  of  refraction  n  is  a  function  of  both  range  r  and  depth,  and 
that  the  method  has  second  order  accuracy  in  r.  Thus,  theoretically 
if  the  range  step  dr  is  sufficiently  small,  the  inherent  error  in  the 
solution  can  be  made  negligible.  Finally,  Jensen  and  Krol  C 193  state 
that  the  PE  method  yields  a  wave  solution  that  includes  all  mode 
coupling  and  diffraction  effects. 

The  PE  method  is  able  to  handle  sloping  bottoms.  By 
partitioning  the  range  into  increments  in  which  the  SVP  is  range- 
independent,  and  then  piecing  these  range  increments  together  via  the 
split-step  or  IP'D  technique,  the  PE  method  is  able  to  numerically 
handle  the  range-dependent  case. 

Diffraction  and  all  other  full-wave  effects,  such  as  the 
rigorous  treatment  of  caustics,  are  included  in  the  PE  method.  Using 
the  standard  PE,  discrete  modes  are  propagated  with  the  correct 
amplitudes  and  mode  shapes,  but  with  errors  in  the  phase  and  group 
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velocities.  These  errors  are  essentially  eliminated  by  the  wide-angle 

PE. 

A  range-dependent  environment  may  be  treated  by  the  PE  method, 
by  a  ray  model,  or  by  the  adiabatic  approximation  to  the  normal  mode 
model.  Of  these,  the  first  two  are  generally  more  accurate  than  the 
third.  The  computation  time  required  increases  with  frequency  for 
wave  models  <mode,  fast  field,  PE),  while  frequency  doesn’t  effect  the 
time  for  ray  models.  The  PE  model  requires  essentially  the  same 
calculation  time  for  the  range-dependent  and  range-independent  cases, 
while  ray  and  mode  models*  time  consumptions  are  proportional  to  the 
number  of  profiles. 

X.  Possible  Disadvantages  of  the  Parabolic  Equation  Method 

Puraoci  vo  rmrotfat  %  on  Hues  reauired  in  a  shallow  water 
environment  have  been  a  major  drawback  for  the  PE  method.  However, 
Jensen  and  Kupernvan  C  23  3  report  increases  in  speed  by  a  factor  of  100 
when  the  PE  model  is  run  on  a  dedicated  computer  system  in  connection 
with  an  array  processor,  as  opposed  to  a  general -purpose  computer. 
This  obviously  makes  the  PE  a  more  practical  alternative  when 
calculating  propagation  loss  in  shallow  aater. 

Computation  time  for  the  PE  method  increases  with  the  square 
of  the  frequency,  so  that  it  is  impractical  for  higher  frequencies. 
Another  disadvantage  is  that  is  difficult  to  handle  shear  propagation 
in  the  bottom.  Also,  the  PE  method  is  thusfar  unable  to  handle 
scattering.  A  possible  method  of  handling  volume  scattering  in 
acoustic  propagation,  by  means  of  an  ordinary  differential  equations 
approach  to  the  PE,  is  offered  by  McDaniel,  Saad  and  Lee  in  a  recent 
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paper  C323.  However,  backscatter ing  must  be  negligible  for 
application  of  the  PE  method. 

Jensen  and  Krol  C  19  3  observe  that  steep  sound-speed  gradients, 


such 

as  those  that 

occur  at  the 

water /bottom 

interface  in  shallow 

water. 

cause  excessive  computer 

run-t imes 

for 

the  PE 

method. 

From 

this. 

together  with 

the  claim 

that  the 

PE 

method 

handles 

range 

dependence  in  a  simple  way,  they  conclude  that  the  primary  application 
of  the  method  should  be  *  low-f requency  acoustic  modeling  in  deep  water 
(no  bottom  effects)  under  range- depen dent  environmental  conditions.* 
Although  the  application  of  the  PE  method  to  shallow  water  propagation 
and  the  incorporation  of  bottom  effects  has  been  improved  since  1975, 
as  described  elsewhere  in  this  paper,  we  question  whether  the  method 
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this  regard  are  described  in  detail  above. 

As  pointed  out  elsewhere  in  this  paper,  the  small  range  steps 
necessary  in  the  PE  method  when  there  is  a  large  vertical  sound  speed 
gradient,  such  as  in  a  shallow  water  environment  at  the  water/bottom 
interface,  result  in  excessive  computation  times.  Also,  as  the  sice 
of  the  range  increment  in  a  numerical  impl imentat ion  of  the  PE  method 
is  inversely  proportional  to  the  frequency,  the  method  rapidly  becomes 
too  time-consuming  as  the  frequency  increases.  At  the  current  state 
of  the  art  in  coding  for  the  PE  method,  the  upper  frequency  limit, 
when  the  program  is  run  on  a  computer  such  as  the  VAX  11/780,  is  about 
500  Hz.  On  a  super -computer  this  upper  limit  could  obviously  be 
raised,  while  on  a  microcomputer  the  highest  practical  frequency  would 
be  much  less  than  500  Hz.  However,  if  an  efficient  dedicated 
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mi croprocessor  could  be  developed,  it  is  conceivable  that  frequencies 
as  high  as  1  kHz  could  be  handled  by  the  PE  method  without  inordinate 
difficulty. 

A  starter  solution,  obtained  from  normal  mode  theory  or  some 
other  method,  is  required  in  order  to  initialize  the  PE  method.  This 
can  significantly  increase  the  computation  time  required.  Also,  the 
reflection  that  can  occur  due  to  a  range-depandent  environment  cannot 
be  accounted  for  by  the  PE  method,  since  all  energy  is  assumed  to 
propagate  in  the  forward  direction. 

The  direct  modeling  of  boundaries  is  not  allowed  for  in  the  PE 
method.  Surface  modeling  is  accomplished  by  making  the  modal  depth 
function  reflected  about  the  surface  with  a  180*  phase  change.  Bottom 
absorption  is  accounted  for  by  modeling  the  bottom  as  a  change  in 
sound  speed  with  an  imaginary  part. 

Ellis  and  Chapman,  in  a  recent  and  very  interesting  paper 
CIA  3,  describe  and  combine  all  extensions  to  the  Pekeris  model,  which 
has  proven  to  be  a  very  useful  tool  in  the  study  of  low-frequency 


propagat ion 

in  shallow 

water.  This 

simple 

model 

contains 

all  the 

basics  of 

norma 1  mode 

theory,  and 

does 

not 

deal  with 

the  PE. 

Nevertheless  their  results  are  of  interest  for  the  current  study,  both 
as  a  means  of  comparison  with  PE  results,  and  for  the  informative 
light  which  they  shed  on  shallow  water,  low-f requency  propagation.  In 
particular,  the  model  accurately  predicts  the  very  high  losses  at 
frequencies  less  than  200  Hz  resulting  from  shear  waves  in  a  chalk 
bottom.  The  reason  that  the  modified  Pekeris  model  performs  much 
better  in  this  environment  than  the  usual  normal  mode  models  appears 
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to  be  that  the  shear  wave  speed  in  the  bottom  is  roughly  of  the  same 
order  of  magnitude  as  the  speed  of  sound  in  mater,  whereas  other 
normal  mode  models  require  that  the  shear  wave  speed  be  either  much 
less  or  much  greater  than  the  water  sound  speed.  Whatever  the  reason, 
the  modified  Pekeris  model  gives  much  better  results,  under  these 
environmental  conditions,  than  any  model  based  upon  the  PE  that  we 
have  seen.  Hence,  it  should  be  considered  as  a  very  desirable 
alternative  to  the  PE  method  for  low-frequency  shallow  water 
propagation. 

XI.  Recommendations  for  Future  Work 

1.  Development  of  a  dedicated  microprocessor,  utilizing  the 
PE  method  to  solve  the  wave  equation. 

2.  Determination  of  the  speed  and  accuracy  of  the  Pesogen 
(Parabolic  Equation  Solution  Generator)  system,  designed  for  shipboard 
use  by  Tappert,  Phu  and  Daubin.  Information  on  this  system  seems  to 
be  unavailable  in  the  open  literature. 

3.  Estimate  the  errors  inherent  in  the  PE,  and  in  the  other 
adaptations  of  the  acoustic  wave  equation. 

4.  Perform  a  sensitivity  analysis,  to  determine  the  effect  of 
errors  in  the  input  parameters  upon  the  solution  of  the  PE. 

3.  Develop  a  method  capable  of  handling  a  rough  surface,  such 
as  will  be  encountered  in  an  underice  environment. 

6.  Develop  a  method  capable  of  handling  scattering,  such  as 
will  be  encountered  in  an  under  ice  environment. 

7.  Develop  a  very  wide  angle  solution. 

8.  Develop  a  broadband  capability. 
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9.  Improve  the  capabilities  of  PE  models  to  handle  three- 
dimensional  problems. 

10.  Incorporate  the  consideration  of  fronts  and  eddies  in  the 
PE  models. 


11.  Develop  methods  of  more  accurately  measuring  the  sound 
speed  profile  in  an  under ice  environment. 

12.  Develop  more  accurate  methods  of  modeling  and  handling 
surface  and  bottom  boundary  conditions. 

13.  Develop  automatic  step-size  determination. 

14.  Develop  the  ability  to  handle  shear  waves. 

15.  Include  multiple  irregular  interfaces  in  a  practical 

model. 


XII.  Conclusions 

With  the  current  state  of  the  art  in  propagation  loss 
modeling,  it  appears  hopeless  to  expect  accurate  absolute  measures  of 
loss  levels  except  in  the  very  simplest  of  environments.  However,  it 
is  quite  likely  that  reasonably  good  comparative  measures  are 
attainable.  For  this  to  be  useful  when  conducting  a  search,  it  is 
necessary  to  have  detailed  knowledge  of  the  environment  in  which  the 
search  is  being  carried  out,  so  that  a  basis  for  a  valid  comparison 
will  exist. 

As  described  in  the  body  of  this  paper,  the  parabolic  equation 
method  will  often  be  the  method  of  choice  when  making  propagation  loss 
predictions.  However,  as  we  have  seen,  a  knowledge  of  the  environment 
is  again  crucial,  since  other  methods  are  preferable  under  those 
circumstances  described  previously. 


£-3A- 


In  general,  good  agreement  between  theory  and  experimental 
data,  whichever  acoustic  model  is  being  employed,  requires  accurate 
knowledge  of  such  environmental  details  as  the  SVP,  bottom  layering, 
bottom  rigidity  (shear),  scattering,  etc.,  as  well  as  the  range- 
dependency  of  these  conditions.  However,  the  required  environmental 
information  is  often  unavailable.  Furthermore,  as  acoustic  models 
become  mare  sophisticated  they  require  even  more  extensive  knowledge 
of  the  propagation  environment.  We  believe  that  the  primary 
difficulty  today  in  ocean  acoustic  modeling  is  the  lack  of 
sufficiently  detailed  and  accurate  environmental  information  for  input 
to  the  models.  To  us  this  appears  more  crucial  than  further 
refinements  to  the  models  themselves.  As  environmental  details  are 
generally  difficult  and  expensive  to  obtain,  a  knowledge  of  the 
accuracy  required  in  the  input  data  to  the  various  models,  in  order  to 
obtain  sufficiently  accurate  propagation  loss  results,  would  seem  to 
be  extremely  important.  An  analysis  of  this  crucial  question  is 
strongly  recommended. 
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Abstract 

A  nearly  ideal  notch  filter,  employing  coefficients  of  equal 
magnitude,  is  described.  Applications  to  the  design  of  transmitting 
antenna  arrays  are  discussed  briefly.  The  construction  is  based  upon 
earlier  work  of  the  author  involving  polynomials  with  restricted 
coefficients. 
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Figure  Captions 

Figure  1.  A  graph  of  for  N=60,  in  dB’s. 
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Abstract 

A  nearly  ideal  notch  filter,  employing  coefficients  of  equal 
magnitude,  is  described.  Applications  to  the  design  of  transmitting 
antenna  arrays  are  discussed  briefly.  The  construction  is  based  upon 
earlier  work  of  the  author  involving  polynomials  with  restricted 
coefficients. 

I.  INTRODUCTION 

The  classical  mathematical  problem  in  notch  filter  design  is 
to  produce  a  polynomial  whose  magnitude  on  the  unit  circle  is  close  to 
constant  in  almost  all  directions,  but  which  has  a  small  number  (i.e., 
I,  2  or  3)  of  deep  nulls  ('notches")  at  specified  points.  Such 
filters  are  applied,  for  example,  to  remove  spectral  lines  from 
otherwise  broadband  spectra.  In  this  paper  we  produce  such 
polynomials  having  one  null,  with  the  added  feature  that  all 
coefficients  have  the  same  magnitude.  For  convenience  this  magnitude 
is  assumed  to  be  one.  Observe  that  this  “unimodular "  property  allows 
the  direct  application  of  these  polynomials  to  the  design  of 
transmitting  antenna  arrays  which  are  omnidirectional  except  for  1 
null.  This  feature  is  crucial  in  certain  communications  areas,  where 
it  is  desired  to  null  out  one  listener  in  a  known  direction  while  at 
the  same  time,  for  maximum  efficiency,  all  antenna  elements  are 
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broadcasting  at  full  power. 

If  the  polynomial  P(z)  is  of  degree  n-1,  it  is  clear  from  the 
Parseval  Theorem  that  its  L2  norm  (i.e.,  RMS  value)  is  exactly  nlx2, 
since  there  are  n  coefficients  each  of  magnitude  1.  Thus,  for  |p(z)l 
to  ba  close  to  constant  on  Jzl  =  l,  that  constant  must  be  n1''2.  The 
question  of  the  existence  of  such  polynomials  is  a  classic  one  in 
mathematical  analysis.  Its  study  uias  apparently  initiated  by  Hardy 
C  11, p.  1993,  and  furthered  by  Littlewood  C8,93,  Erdos  C5  3,  Newman 
C 1,2,3, 10]  and  others.  A  basic  result  concerning  these  problems  was 
obtained  by  the  author  [43,  which  paved  the  way  for  solutions,  by 
Korn»r  C73  and  Kahane  C6  3,  of  two  of  the  fundamental  conjectures  in 
this  area.  We  modify  the  construction  given  in  [43  to  produce  nearly 
ideal  filters  with  one  notch. 


by 


II.  APPROACH  AND  RESULTS 

Our  starting  point  is  the  polynomial  P,  of  degree  N*-l,  given 


p<e>= 

M-l  M'-/ 

Y  /  exp  (2TTi  jkN“M  z=*exp  (2rTiO). 

It  is  shown  in 

t=0  j--0 

[43  that  P(-O)  satisfies! 

( i  ) 

|  P  ( roN~a  >|=N  for  all  integers  m, 

<lie/<i/2. 

( ii ) 

For  any  £,  N“‘<«l/2,  Jp(0)l=N+E  for 

{  iii) 

IeI  <l+2TT-*+5<Tr€)-S 

For  N  odd,  P< 1/2N)=0< 1 > ,  while 

for  N  even 

1 >/<2N2) )=0< 1 ) 

,  and 

where 


v  P( <N- 
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(iv)  lp<«)l<(2+3rr-\>N  +  0(1)  for  all  e. 

Recent  numerical  evidence  suggests  that  (iv)  can  be  strengthened  to: 
(iv’)  jp(-0)l<1.4N  for  all  -0. 

Thus,  as  N-»oo,  the  magnitude  of  P  is  asymptotically  close  to  the  ideal 
constant  N  except  for  the  immediate  neighborhood  of  one  point.  By  a 
simple  change  of  variables  it  is  clear  that  this  special  point  can  be 
taken  anysuhere  on  the  unit  circle. 

If  P(O)  is  changed  by  removing  the  first  N  terms,  all  of  whose 
coefficients  are  +1,  and  then  dividing  by  zM,  there  results 


A/-/  hh  I 


Q(-e>= 


L  L 


e x p  ( 2TT  i  j  kN~ 1 )  z  «• ♦ 


This  Q,  which  is  of  degree  N2-N-l,  is  the  desired  modification  of  P. 
In  fact,  for  4>=0  the  terms  in  the  inner  sum  for  Q  are  identically  0, 
so  that  Q ( 0 >  =0 .  Also,  the  N  terms  removed  from  P  are  significant  only 
in  the  immediate  neighborhood  of  0=0  (i.e. ,  the  special  point),  so 
that  estimates  (ii)  and  (iv’)  remain  true  for  Q.  Furthermore,  it  can 
be  shown  that  the  null  width  of  Q  is  less  than  2/N.  These  properties 
together  with  figure  1,  which  exhibits  1 Q (O > i  in  dB’s  as  a  function  of 
•O  for  N=6@  (Q  of  degree  3539),  show  that  Q  is  indeed  the  nearly  ideal 
notch  filter  discussed  earlier.  Once  again  it  is  clear  that  a  change 
of  variables  allows  the  relocation  of  the  notch  to  any  desired  <K 
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ABSTRACT 

Several  new  designs  of  analytic  null  steering  algorithms  for 
linear  arrays  are  described.  Two  of  them,  the  ^-Technique  and  the 
Positive  Coefficient  Model,  allow  for  placing  an  arbitrary  number  of 
nulls  in  arbitrary  directions,  while  maintaining  main  beam  and 
Sidelobe  level  control.  A  method  of  incorporat ing  these  deterministic 
null  steering  techniques  into  existing  adaptive  algorithms,  offering 
the  possibility  of  increases  in  array  performance  at  small  cost,  is 
currently  being  researched. 
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Nall  Steering  Employing  Polynomials  with  Restricted  Coefficients 
J.  S.  Byrnes  (Sr.  Member,  IEEE)  and  Donald  J.  Newman 

Prometheus  Inc. 

I.  INTRODUCTION 

In  view  of  the  well-known  one-to-one  correspondence  between 
polynomials  and  linear  arrays  with  commensurable  separations  between 
elements,  as  described  in  detail  by  Schelkunoff  C93,  null  steering 
questions  involving  such  arrays  translate  directly  into  mathematical 
problems  regarding  the  locations  of  zeroes,  on  the  unit  circle,  of 
polynomials*  Furthermore,  physical  and  electronic  limitations  placed 
upon  the  array  elements,  such  as  a  maximum  allowable  transmitting  power 
or  a  bound  on  the  dynamic  range,  imply  various  restrictions  upon  the 
coefficients  of  these  polynomials.  Here,  dynamic  range  refers  to  the 
ratio  of  the  magnitudes  of  the  largest  to  the  smallest  weight,  or 
Shading  coefficient,  of  the  array.  Thus,  the  theoretically  challenging 
question  of  the  placement  of  zeroes  at  specified  points  on  the  unit 
circle,  of  polynomials  whose  coefficients  satisfy  certain  restrictions, 
is  also  a  problem  of  strong  practical  interest  to  antenna  designers. 

The  design  af  filters  is  another  application  in  which  such 
questions  arise.  For  example,  the  classical  mathematical  problem  in 
notch  filter  design  is  to  produce  a  polynomial  whose  magnitude  on  the 
unit  circle  is  close  to  constant  in  almost  all  directions,  but  which  has 
a  small  number  (ie,  1,  2  or  3)  of  deep  nulls  (’notches')  at  specified 
paints.  In  C53  the  construction  of  C43  is  employed  to  produce  such  a 
polynomial  having  one  null,  with  the  added  feature  that 
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all  coefficients  have  the  same  magnitude.  Hence,  the  dynamic  range  of 
the  notch  filters  presented  in  £5  3  is  one. 

This  paper  addresses  the  null  steering  application  described 
in  the  first  paragraph.  There  are  several  factors  which  must  be 
considered  in  the  design  of  null  steering  algorithms.  In  addition  to 
the  basic  problem  of  placing  the  nulls  the  main  beam  must  be  steered, 
the  width  of  the  main  lobe  controlled,  and  the  sidelobe  levels  must  be 
sufficiently  below  that  of  the  main  lobe.  Control  of  the  sidelobe 
level  is  usually  achieved  by  attenuating  the  shading  coefficients  as 
one  moves  away  from  the  center  of  the  array.  Often  these  attenuation 
factors  (Chebyshev,  Taylor,  etc.)  are  chosen  in  advance,  and  may  not 
be  easily  altered  once  the  array  is  in  place.  This  leads  directly  to 
a  beautiful  mathematical  question,  similar  to  the  peak-factor  problem 
in  engineering  attacked  by  Boyd  £23,  Schroeder  £103  and  others: 

Given  the  magnitude  of  the  coefficients  of  a  polynomial  P,  a 
finite  subset  3  of  the  unit  circle  C,  and  a  point  p<C  distinct 
from  those  in  S,  choose  the  phases  of  these  coefficients  so 
that  P(z>=0  for  all  z<S,  the  maximum  on  C  of  1 P ( z ) 1  occurs  at 
z=p,  and  the  maximum  of  Jp<z)J  on  a  subset  of  C  excluding  an 
appropriate  interval  (the  beamwidth)  around  p  is  as  small  as 
possible. 

We  consider  various  subprablems  in  this  paper.  Research  on  the 
general  question  is  continuing. 

II.  DIRECT  ADAPTIVE  NULLING  1  ' 

Currently  the  most  widely  used  class  of  null  steering  methods 
is  known  as  adaptive  nulling  £1,3,6,7,8,123.  Adaptive  arrays  have 
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developed  over  the  past  tuienty-five  years  as  the  preferred  method  of 
reducing  the  performance  deterioration  in  signal  reception  systems  which 
is  inevitably  caused  by  noise  entering  the  system.  These  methods,  such 
as  the  Widrow  LMS  Algorithm  and  the  Howel ls-Applebaum  sidelobe 
canceller,  are  indirect  adaptive  schemes,  they  do  not  explicitly  form  an 
estimate  of  the  directions  of  arrival  of  interfering  sources  or 
explicitly  steer  nulls  in  those  directions.  A  scheme  in  which  these  two 
tasks  are  actually  performed  can  be  called  a  direct  adaptive  algorithm. 
Thus  one  possible  approach  to  the  solution  of  many  of  the  well-known 
performance  problems  which  arise  with  existing  methods  is  to  complement 
an  appropriate  indirect  adaptive  algorithm  with  the  analytic  null 
steering  methods  described  herein.  In  this  way  the  actual  noise 
suppression  achieved  might  be  enhanced  beyond  that  which  would  be 
available  through  either  adaptive  or  analytic  methods  exclusively.  In 
addition,  since  each  execution  of  our  analytic  algorithms  is  essentially 
instantaneous,  and  since  a  direct  scheme  allows  much  greater  use  of 
prior  knowledge,  such  as  known  jammer  locations  or  known  multipaths, 
this  interactive  direct  method  offers  the  possibility  of  increases  in 
performance  at  small  cost.  Research  on  this  approach,  which  we  call  the 
’'Direct  Adaptive  Nulling  System,*  is  in  progress. 

III.  COEFFICIENTS  OF  EQUAL  MAGNITUDE 
An  important  subproblem  of  the  general  mathematical  question 
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described  earlier  is  the  case  when  all  of  the  coefficients  of  the 
polynomial  have  the  same  magnitude,  which,  by  normalization,  we  can 
assume  to  be  one.  Such  phase-only  shading  occurs,  for  example,  in  the 
design  of  transmitting  arrays  which  are  omnidirectional  except  for 
specified  nulls.  These  features  are  crucial  in  certain  communications 
areas,  where  it  is  desired  to  null  out  listeners  in  known  directions 
while,  at  the  same  time,  for  maximum  efficiency,  all  antenna  elements 
are  broadcasting  at  full  power.  Also,  in  order  to  minimize  the 
relative  size  of  the  quantisation  steps  in  a  gradient  algorithm  such 
as  LMS,  the  coefficient  magnitudes  should  be  kept  as  close  as  possible 
to  unity  L 7,  p.  1531.  Note  that  this  “equimagnitude"  property  of  the 
coefficients  precludes  the  use  of  attenuators,  with  the  concomitant 
savings  in  electronic  hardware. 

The  most  elementary  example  of  the  above  is  the  unshaded  array 
-  all  coefficients  are  1.  In  spite  of  its  simplicity,  this  uniform 
array  is  of  practical  importance.  Observe  that  in  this  case  the 
zeroes  of  the  polynomial  are  almost  uniformly  spaced  around  the  unit 
circle,  occurring  at  all  of  the  n+1**  (where  n  is  the  degree  of  the 
polynomial)  roots  of  unity  except  2=1,  where  there  is  a  maximum. 

At  the  rather  extreme  is  the  case  where  an  n-fold  zero  _is 
required  at  one  point.  One  application  of  this,  as  discussed  by 
Steyskal  till,  is  to  broaden  a  pattern  null  so  as  to  null  an  entire 
sector.  Clearly,  by  a  simple  change  of  variables,  this  zero-point  can 
be  assumed  to  occur  at  z=l.  It  is  a  straightf crwarcl  matter  to 
construct  such  a  polynomial  with  coefficients  of  magnitude  1;  in  fact, 
the  coefficients  may  all  be  taken  to  be  +1.  Namely,  define  P(z)  by: 


n-1 

P(  z)='n(l-zs*m)  (1) 

m=0 

The  problem  with  this  construction  is  that,  although  PCz)  obviously 
satisfies  the  required  properties,  it  does  so  at  very  high  cost. 
Since  P  has  degree  2n-l ,  its  realization  requires  an  array  with  2n 
elements.  We  show  in  the  following  theorem  that,  for  all  but  small 
values  of  n,  this  situation  may  be  greatly  improved  by  allowing 
coefficients  to  be  0  as  well  as  +1.  Since  this  simply  means  that  some 
array  elements  are  turned  off,  the  dynamic  range  of  the  coefficients 
is  not  affected  in  any  meaningful  way. 

Theorem  l  Let  n>10.  Then  there  is  a  polynomial  P(z),  of 
degree  less  than  n3,  such  that  P(z)  has  an  n-fold  zero  at  z=l,  and  all 
the  coefficients  of  P  are  either  +1  or  0. 

Proof  of  Theorem  J[  Given  n>10,  choose  k  so  that 

2fcl>|<”2  (2). 

Since  x/ln(x)  is  increasing  for  x>e,  and  lE/ln ( 10)>3/ln(2) ,  the  choice 


k=n3  certainly  implies 

(2) . 

Actually,  for 

any  O0,  it  is  clear 

that 

for  n  large  enough  we 

may 

take  k^n3"^  ,  but 

such  precision 

is 

not 

necessary  here.  Now, 

for 

each  arbitrary 

subset  S  of  the 

set 

of 

nonnegative  integers  less  than  k,  let 


Q(z>  =  ZL  z~, 
»uS 
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and  form  the  vector 


(Q(l)t  Q’(l),  Q,Ml)/2!,...Q,"-‘,(n/tn"l)!). 

These  are  integer  vectors*  and  the  largest  entry  is  bounded  by 

k-i 

Y  k". 

rr\zO 

t 

Thus  there  are  less  than  k"  such  vectors.  Since  there  are  2**  subsets 
S,  and  so  2 “  polynomials  Q,  equation  (2)  implies  that  at  least  two 
distinct  polynomials,  say  Qi(z)  and  Qa(z),  have  the  same  associated 
vector.  Hence  PtzJ^CUtzJ-Qatz)  is  the  desired  polynomial,  and  Theorem 
I  is  proven. 

The  idea  underlying  equation  (1),  which  we  call 
“encapsulation, “  may  also  be  employed  to  construct  polynomials  with 
coefficients  of  magnitude  1  that  place  any  number  of  arbitrary  nulls. 
Namely,  we  have: 

n 

Ijjegrgm  II  For  any  positive  integer  n,  let  Cz«3  be  an 

*^•1 

arbitrary  set  of  (not  necessarily  distinct)  points  on  the  unit  circle. 
Then  there  is  a  polynomial  P(z>  with  coefficients  all  of  magnitude. 1, 
of  degree  2*'~1,  satisfying  P(zm)=0,  l<m^n. 

Prqqf  of  Theorem  IX  As  indicated  above,  we  simply  produce  an 
explicit  formula  for  P(z): 

*  v 

t\  *W  h<| 

PCz)  =  TT  (za  -z«a  ). 
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It  is  straightforward  to  see  that  this  P(z>  satisfies  the  required 
properties,  establishing  Theorem  II. 

IV.  SIMULTANEOUS  NULL  STEERING  AND  MAIN  BEAM  PLACEMENT 

In  the  previous  section  we  attacked  the  subproblem,  of  the 

general  question  stated  earlier,  which  arises  when  the  dynamic  range 

of  the  coefficients  is  required  to  be  one.  We  now  consider  another 

aspect  of  the  original  problem,  perhaps  of  more  interest  to  antenna 

designers.  Namely,  how  can  arbitrary  nulls  be  placed  while 

maintaining  a  specified  main  beam  direction  and  specified  maximum 

sidelobe  level?  We  describe  two  methods,  the  ^-Technique  and  the 

Positive  Coefficient  Model,  of  achieving  these  goals. 

To  set  the  problem,  again  let  n  denote  any  positive  integer, 
n 

and  let  S«Cz^3  be  an  arbitrary  sst  of  (not  necessarily  distinct) 
points  on  the  unit  circle.  Also,  let  z*  be  a  point  on  the  unit  circle 
distinct  from  those  in  S.  Our  methods  allow  the  placement  of  zeroes 
of  a  p  ynomial  P  at  all  points  in  S,  while  simultaneously  having  the 


maximum  of 

1  Pi  on 

the  unit 

circle  occur  at 

z=z».  Furthermore, 

-  the 

difference 

between 

iP(z«)i 

and 

the  highest 

sidelobe  can  be 

made 

arbitrari ly 

large. 

As  will 

be 

seen,  the 

costs  encountered 

in 

achieving  the  last  property  are  an  increase  in  the  degree  of  P,  and  a 

loss  of  control  of  the  dynamic  range  of  the  coefficients. 

To  proceed  with  the  constructions,  define  the  angles  £0«3, 

by  z««®exp(  ,  -Ik0«^ir.  As  before,  a  simple  change  of 

variables  allows  us  to  assume  O«=0,  so  that  z«®l. 

Method  1.  The  ^-Technique 
n 

Let  i=-cot~*C2  cot<^«,/2)  3,  z*=exp<  iS) , 

/ 
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and  define  Q(z>  by 


n 

Q(z)=(z-z*)TT  ( z-zm) . 

"\il 

A  straightforward  calculation  shows  that  1q(z)I  has  a  relative  maximum 
at  z=l .  Hence,  for  c  a  large  enough  positive  integer, 
PtzJsd+zJ^QCz)  will  certainly  satisfy  the  required  properties.  It 
can  be  shown  that,  in  order  to  guarantee  an  absolute  maximum  of  lp(z)J 
at  z=l,  it  is  sufficient  to  take  c*=(i/<],  where  <=min/o,*/.  Of  course, 
in  order  to  further  increase  the  main  lobe  level  relative  to  the 
sidelobes,  if  will  be  necessary  to  take  c  larger. 

Method  II.  The  Positive  Coefficient  Model 

For  each  m,  l<m<n,  choose  the  smallest  positive  integer  k* 
such  that  exp(ik„,<5u)  lies  in  the  left  half  plane,  and  define  P(z)  by 


Clearly  P(z 

)  has 

the 

necessary 

zeroes. 

Furthermore, 

al  1  of 

the 

coefficients 

of  P 

are  po 

sitive, 

so 

that  the 

maximum  of 

1p(z)1  on 

the 

unit  circle 

obviously 

occurs 

at 

z=»l. 

Once  again  it 

is  simple 

to 

further  increase 

the  main  lobe 

level  relative  to  the 

sidelobes, 

by 

multiplying  P(z)  by  an  appropriate  positive  integer  power  of  l+z. 

There  are  two  additional  points  which  can  be  made  about  the 

.  •  * 

Positive  Coefficient  Model.  One  is  that  its  electronic  implementation 
will  be  greatly  simplified  as  compared  to  that  of  arbitrary  shading 
coefficients,  since  the  positivity  eliminates  the  need  for  phase 
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shifters.  A 

second 

is 

that  some 

control  of 

the 

dynamic  range  of 

the 

coefficients 

can  be 

achieved  by 

combining 

this 

method  with 

the 

encapsulation 

technique 

discussed 

earlier, 

if 

we 

again  ignore 

the 

effects  of  0  coefficients. 

V.  CONCLUSION 

Various  aspects  of  a  fascinating  problem  in  classical 
mathematical  analysis,  with  direct  applications  to  antenna  array  design, 
have  been  discussed,  and  several  results  obtained.  Foremost  among  these 
are  two  analytic  methods  for  placing  an  arbitrary  number  of  nulls  in 
arbitrary  directions,  while  maintaining  main  beam  and  sidelobe  level 
control.  A  method  of  incorporating  these  analytic  null  steering 
techniques  into  existing  adaptive  algorithms,  offering  the  possibility 
of  increases  in  array  performance  at  small  cost,  is  currently  being 
researched. 
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Abstract 

Precise  estimates  were  previously  obtained  for  the  error  in 
the  optimum  gain  of  an  antenna  array  arising  from  errors  in  the 
measured  noise  covariance  matrix.  It  was  observed  that  the  error 
estimates  are  a  function  of  the  character  and  spatial  distribution  of 
the  noise  field.  The  purpose  of  the  current  work  is  to  analyze  this 
dependence  of  the  previous  results  upon  the  quantitative  and 
qualitative  aspects  of  the  noise  field. 


A.  Introduction 


In  C13  the  question  of  the  effect  of  errors  in  the  measured 
noise  covariance  matrix  upon  the  gain  of  an  antenna  array  mas  addressed. 
Precis*  error  estimates  were  obtained  for  the  case  of  a  linear  array  of 
n  equally  spaced  elements  in  a  noise  field  made  up  of  spatially 
uncorrelated  point  sources  of  a  single  frequency.  Note  that  under  these 
conditions  the  noise  covariance  matrix  is  Hermitian,  positive  definite 
and  Toeplitz.  It  was  observed  in  Ill  that  the  error  estimates  were  a 
function  of  the  character  and  spatial  distribution  of  the  noise  field. 
The  purpose  of  the  current  work  is. to  analyze  this  dependence  of  the 
previous  results  upon  the  quantitative  and  qualitative  aspects  of  the 
noise  field.  Several  of  the  results  given  below  ware  announced  in  C23. 
A  significant  improvement  in  the  major  result  (Theorem  C),  as  well  as 
proofs  for  the  previously  announced  theorems  and  a  discussion  of  their 
importance,  is  given  herein. 


B.  Results 

The  notation  will  be  the  same  as  that  in  C 1 3  and  C2  3,  so  that  n 

is  the  number  of  antenna  elements,  R  is  the  true  noise  covariance 

a  a 

matrix,  R  is  the  measured  noise  covariance  matrix,  E^R-R  is  the  error 

»atrix,ll*ti  is  a  norm  on  the  space  of  all  n  x  n  matrices,  and  i«*lisU. 

Far  weight  <i,e.  "shading")  vector  W  and  signal  vector  S  the  gain  is 

given  by 


G* 


|(Wf,S)j  » 
- 

U  RU 


When  U  is  chosen  to  be  an  ideal  Wiener  weight  vector, 


U»R— ‘Sv 


the  resulting  6  is  known  as  the  optimal  gain.  Since  the  measured  noise 

A 

covariance  matrix  R  is,  by  necessity,  employed  in  obtaining  the  actual 
Weiner  weight  vector 

a  a 

W*R—‘S, 

the  actual  gain  will  differ  from  the  optimal  gain  G.  Denoting  this 

actual  gain  by  G,  it  is  seen  that 

t  A 

(3  R~*S)a 

r— 1 * 

t  ^  A 

S  R"1RR“tS 

Note'  that  the  true  noise  covariance  matrix  R  appears  in  the  denominator 
of  6. 

In  order  to  determine  the  effects  of  the  noise  field  upon  the 
gain,  the  following  result  from  C 1 3  will  be  employed: 

Theorem  A  Cl,  Corollary  III.  For  any  small  enough  positive  number 

G 

—  -  i  <  16||r-‘||*||rH*. 

G 

It  is  clear  from  Theorem  A  that  what  is  now  required  is  an 
analysis  of  how  the  norms  of  the  noise  covariance  matrix  and  its  inverse 
depend  upon  the  noise  field.  As  is  well  known,  the  determination  of  the 
norm  of  a  matrix,  except  in  very  special  cases,  is  no  easy  mattor.  Even 
for  the  particular  class  of  matrices  considered  here  (i.e.,  Hermitian, 
positive  definite,  Toeplitz),  precise  estimates  for  the  norm  are 
difficult  to  attain.  The  most  tractable  norm  to  deal  with  in  this 


context  is  defined,  for  an  n  x  n  raatrii;  A,  by: 


IHI*  -  Max  jy~Axj , 


where  the  maximum  is  taken  over  all  n-vectors  x  and  y  with  Euclidean 
norm  1.  If  A  is  Hermitian  and  positive  definite,  and  if  ,Xt»  ■  ■  •  »^,i» 
are  the  (necessarily  real  and  positive)  eigenvalues  of  A  written  in 
increasing  order,  it  follows  from  the  standard  theory  that 

i|Aj]a  =  X  and  .  .  . 

n 

<  i) 

I|a— 


Hence,  accurate  estimates  for  the  smallest  and  largest  eigenvalues  of  R 
ar®  required. 

Consider  a  model  of  an  isotropic  noise  field  plus  a  single  point 
source  of  noise  which  is  generating  a  plane  wave  signal.  Let *  be  the 
power  of  the  isotropic  noise,  g^  the  point  source  power,  \  the 
wavelength  of  the  noise,  and  assume  the  array  spacing  to  be  /C /2.  In 
addition,  define  the  ratio  7*  by: 


<T-  °;v 


(2) 


As  is  well  known  CLeBlanc,  p.  12,  for  example 3,  the  normalised  noise 
covariance  matrix  under  these  conditions  is  given  by: 


I 


K 


I 


r 

w 


R» 


/  t  r 
r  i  T 


r  r  t 


•  r 

•  r 


•  « 


A  straightforward  computation  shows  that 


Eigenvalues  ( R)»l-T,  1-Y,  . . . , l~Ti  l+<n— 1  >T, 
so  that  llRll-lXn-1  )T  and  ■Ir~*//=  _ ! - 

11  »-T 


If  the  model  is  altered  by  changing  the  wavelength  of  the  point 
source,  so  that  the  array  spacing  remains  at  one-half  wavelength  of  the 
isotropic  component,  but  not  of  the  point  source,  the  normalized  noise 
covariance  matrix  becomes: 


Tc,  TS. 

'  Tc. 


i 


(5) 


I 


Tc  ...  / 

U.  •**1 

where  c^#‘«“*  *,  d»array  spacing,  tsangle  of  arrival  of  plane  wave 

(t*0  is  endfire),  and  K  is  the  wavenumber  of  the  point  source.  Note 
that  c—mi*cm  *nd  also  that  c„*c*,  where  c^e1*- 

It  is  surprising  that  the  eigenvalues  of  R  are  identical  to  the 
previous  case.  Namely, 

Theorem  B.  Let  R  be  given  by  (3).  Then  the  eigenvalues  and 
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norms  of  R  arm  thosm  appearing  in  <4). 

Proof  of  Theorem  B.  Let  I  denote  the  identity  matrix,  'J  be  the 
matrix  all  of  whose  entries  are  1,  and  V  be  the  matrix  R  of  <3)  with  all 
of  the  X*  replaced  by  1.  Thus  R»< 1-y) I+TV.  In  addition,  V  may  be 
factored,  giving 


10..  *0 

n  0  .  .  .  0 

0  c  •  •  *  0 

'  0  c  •  •  •  0 

i  • 

a 

U  , 

. 

• 

t 

*  Mi 

i 

• 

*•/ 

L  0  0  .  -  ‘  z  l 

0  0  *  •  •  c 

From  this  factor izat ion  it  follows  that 

eigenvalues  (7V)*X  eigenvalues  <U)»0,0, . . . ,0,Tn* 

But  the  above  expression  for  R  implies  that  \  is  an  eigenvalue  of  R  if 
and  only  if  T”l4vL>  is  an  eigenvalue  of  "JV,  and  Theorem  B  is  proven. 

Combining  Theorems  A  and  B  immediately  yields  the  desired  error 
estimate: 

Theorem  C.  Let  R  be  given  by  (3),  and  let  G,  G,  £  ,  X  be  defined 
as  above.  Then,  as  long  as  €.  is  small  enough, 

1 T  -  'I  ‘  «  (TTf  b MT)’ v , 

and  this  estimate  is  best  possible. 

How  consider  a  model  of  two  point  sources  of  noise,  generating 
plane  wave  signals  of  power  and  respectively,  plus  an  isotropic 

noise  field  of  power  Let  KT  and  Kt  be  the  wave  numbers  of  the 

point  sources,  t,  and  tt  the  angles  of  arrival  of  the  plane  waves,  d  the 
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array  spacing,  and  assume  again  that  d^/2,  where  \  is  the  wavelength  of 
the  isotropic  noise.  Defining  T,  d^,  ,  and  f^  by: 

itr\Ktc{  cos  t( 


-r-  Sl±£  -  J  - 
^  +\  f 


‘  M  /T,  c/  Otf  t. 


*<*  -t5—;  >  an^  ^  =  ‘<‘Ci+6-fjeM  , 

<5;  f  (T 


the  normalized  noise  covariance  matrix  is: 


1 

rf, 

TX  *  1 

1  •  Tf„.( 

Tf 

-t 

i 

Tf,  •  • 

-  Tf 

h-  ? 

R* 

Tf, 

• 

9 

Tf 

-( 

/  * 

1  w 

*  •  'Tf  , 

.  *»-3 
• 

Tf 

->  -hrj 

• 

i 

•  » 

• 

1 

(6) 


Theorem  B  can  be  applied  directly  to  obtain  an  upper  bound  for 
1 1  r[  j .  To  see  this,  let  A  denote  the  matrix  R  of  (6)  with  f-,  replaced 


by  c^,  and  let  B  be  R  with  replaced  by  Thus  R^A+d-OB,  so  that 

UrU«*<1Ia[I+  (l-oll»IUn-(n-l)T.  (7) 


Clearly  estimate  (7)  cannot  be  improved  without  additional  assumptions 

t  \ 

regarding  the  relative  sizes  of  (oj  and  0^  ,  so  that  in  the  general  case 
(7)  is  best  passible. 

As  often  occurs  in  estimation  problems,  the  lower  bound  on  the 
eigenvalues  of  R  is  considerably  more  difficult  to  obtain.  It  appears 
necessary  to  employ  asymptotic  estimates,  such  as  the  following  theorem 
CGrenander  &  Szego,  p. 641: 

Theorem  D.  Let  f<x)  be  a  real-valued,  27T-per iodic,  Lebesgue 

am 

intagrable  function  with  Fourier  series  f ( x  )— Zf «exp( imx ) .  (Note  that 
t»"f-.»).  For  each  fixed  n,  n»l,2,3, . . . ,  let  the  (necessarily  Hermitian, 
Toeplitz)  matrix  (f«_«.>,  m, k»l , 2, . . . , n,  have  (necessarily  real) 
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eigenvalues  ^ j * ^* • * -X*  ^ *  Finally,  let  M*  be  the  essential  inf 

of  f<x),  and  let  Ma  be  its  essential  sup  d1t=-<»  and/or  Mas«o  are 
allowed).  Then 


Mi<  yL  <Mav  and 


<8 ) 


lim  lim  =Ma. 

n-*** 


To  apply  Theorem  D  to  estimate  the  smallest  eigenvalue  of  the 
matrix  R  given  in  (6),  first  approximate  R  by  multiplying  each  f^  by  the 
•convergence  factor*  where  0<r<l  and  r  i3  close  to  1,  so  that  now 

f  --Tr1”’ 


>n 


wf]  lff\Kxcastx 

f  O-«0e. 


and 


t 

f&}»  l*~Tc{  l 


m:-40 

!V\f0 


n\:-‘ 
m  £0 


N 


> 

'  I  ^  r  /toj 

'y':’oB  '  '  Ini  /  1  ^  1 

+T(i-4  <L  r  l(^  r^(KxJ“s\+%)+LSnt*(l{<l<ost+'x)l , 

•  2  \  f 

*\io 


Applying  the  standard  formula  for  the  Poisson  kernel,  along  with  the 
facts  that  the  cosine  function  is  even  and  the  sine  function  is  odd,  to 
(9)  yields;  x 

f(x)  =  /+  iTc<[t 


cos  (K,°(  rju^  txj  +  rx  1 

rtTHr^  t  t  -  -1 

yL‘L  l-Zr<os  (Kxc(  cortx  +  x)+  ^ 


(10) 


Clearly,  for  any  values  of  KptMK,t  and  t%,  there  will  be  values  of  x 
such  that  the  first  terms  in  each  of  the  square  brackets  in  (10)  can  be 
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arbitrarily  close  to  0  simultaneously  (for  r  close  enoi-gh  to  1),  so  that 
the  best  possible  lower  bound  for  the  smallest  eigenvalue  of  R  given  by 
(6)  is  1-T.  Applying  (1),  (7)  and  Theorem  A  to  this  yields: 


Theorem  E.  Let  R  be  given  by  (6),  and  let  S,  S,C  ,T  be  defined 
as  above.  Then,  as  long  as  £  is  small  enough. 


-  I  <16 


l 


(,-T)' 
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and  this  estimate  is  best  possible*. 

Next  consider  a  model  where  R  is  a  real,  tridiagonal  matrix. 
Thus,  for  some  real  number  b: 
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As  is  well  known  Cliorgera  and  Cooper,  p.  732  for  example!,  the 
eigenvalues  of  this  R  are: 
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(12) 


Assuming  2b<l,  in  order  to  guarantee  that  R  is  positive  definite,  it  is 
clear  that  the  best  possible  bounds  for  the  eigenvalues  given  in  (12) 
are: 

1  -2b<^L f  < 1 +2b .  (13) 

Combining  (1)  and  (13)  with  Theorem  A  proves: 


V  -a- 


Theorem  F.  Let  R  be  given  by  (11),  assume  2b<l,  and  let  G,  G^C, 
be  defined  as  above.  Then,  as  long  as  £  is  small  enough. 
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and  this  estimate  is  best  possible. 

Finally,  the  complete  eigenspectrum  for  real,  symmetric  matrices 
of  orders  3  and  4,  and  the  partial  eigenspectrum  for  such  matrices  of 
order  5,  have  been  computed  explicitly.  By  various  complex  algebraic 


manipulations,  with 
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the  following  table  may  be  derived: 
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Table  1.  Eigenvalues  and  Eigenvectors  of  Real,  Symmetric  Matrices 


Concentrating  for  the  moment  on  n=3,  assume  that  there  is  a 
real-valued  point  source  of  noise  of  the  form  C/1  cos (Kd  cos  t)  in  an 
isotropic  noise  field  of  power  Defining  "T  by  (2)  and  setting  w=Kd 
cos  t,  the  normal iced  noise  covariance  matrix  is  nous 


I  W  W 
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tW  Joar  yj 


(14) 


The  eigenvalues  of  Rj  may  tie  computed  exactly,  and  the  result  is 
that  they  are  1-]T,  1-yeos  2ui,  and  1+y+ycos  2w.  Applying  (1)  and  Theorem 
A  yieldss 

Theorem  G.  Let  R  be  given  by  (14),  and  let  G,  G,£  ,y  be  defined 
as  above.  Then,  as  long  as  £.  is  small  enough, 


/  <16 


and  this  estimate  is  best  possible  to  within  a  factor  of  4. 

A  similar  result  holds  for  n**4.  The  complete  eigenspectrum  for  n«5  was 
not  derived  in  a  meaningful  form. 

C.  Conclusion 

It  is  clear  from  the  cases  examined  in  Theorems  C  through  G  that 

♦ 

the  angle  of  arrival  of  a  plane  wave  from  a  point  source  of  noise  is  not 
the  crucial  issue,  when  considering  the  accuracy  required  in  measuring 
the  noise  covariance  matrix  of  a  linear  array  of  equally  spaced 
elements.  The  overiding  concern  is  the  relationship  between  the  power 
of  the  point  source  and  the  power  of  the  isotropic  background  noise.  If 


:  -i<3- 


the  ratio  of  the  point  source  power  to  the  isotropic  pokier  is  very  large, 
then  extremely  accurate  measurements  of  the  noise  covariance  matrix  are 
required  in  order  to  get  a  reasonable  estimate  of  the  gain  of  the  array. 
The  same  applies  when  there  are  two  localized  noise  sources  in  an 
isotropic  background,  such  as  occurs  with  torpedo  flow  noise  and 
structural  vibrations.  Nate  the  concurrence  of  these  conclusions  with 
those  of  Hudson,  who  observes  that  the  noise  covariance  matrix  will 
usually  be  ill-conditioned  when  there  is  a  relatively  weak  uncorrelated 
noise  component  in  interference  that  is  due  to  a  small  number  of  discrete 
sources  (less  than  the  number  of  array  elements)  C4,  p. 1361. 
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